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correlations in the cosmological large-scale structure 
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ABSTRACT 

I review the theory of angular momentum acquisition of galaxies by tidal torquing, the result- 
ing angular momentum distribution, the angular momentum correlation function and discuss 
the implications of angular momentum alignments on weak lensing measurements: Starting 
from linear models for tidal torquing I summarise perturbative approaches and the results 
from n-body simulations of cosmic structure formation. Then I continue to discuss the valid- 
ity of decompositions of the tidal shear and inertia fields, the effects of angular momentum 
biasing, the applicability of parameterised angular momentum correlation models and the 
consequences of angular momentum correlations for shape alignments. I compile the result of 
observations of shape alignments in recent galaxy surveys as well as in «-body simulations. 
Finally, I review the contamination of weak lensing surveys by spin-induced shape alignments 
and methods for suppressing this contamination. 

Key words: cosmology: large-scale structure, gravitational lensing 



1 INTRODUCTION 

This article reviews models of angular momenta of galaxies and 
angular momentum correlations in the cosmological large-scale 
structure, and the implications of angular momentum induced 
shape alignments of galaxies for gravitational lensing. The current 
paradigm states that galaxies acquire t heir angular momentum by 
a mechanism known as tidal torquing ([Hovlelll949l ; |Peebleill969l : 
ISciamdl 19551 ; lDoroshkevidJll970l : Iwhitelll984ir where tidal fields 
from the ambient large-scale structure excert a torquing moment 
onto the protogalactic object prior to gravitational collapse. The an- 
gular momentum, which is conserved during collapse, has impor- 
tant implications for the formation and orientation of the galactic 
disk. Tidal torquing explains naturally angular momentum correla- 
tions, as neighbouring galaxies form from correlated initial condi- 
tions and experience similar tidal fields. An important application 
of angular momentum correlations is gravitational lensing: The an- 
gular momentum-induced ellipticity alignments are a source of sys- 
tematics, which, if uncorrected, deteriorates constraints on cosmo- 
logical parameters and introduces estimation biases. 

After a compilation of the key formulas of cosmology, struc- 
ture formation and perturbation theory in Sect. |2l I review angu- 
lar momentum build-up by tidal torquing, which is a second order 
perturbative process arising in Lagrangian perturbation theory, in 
Sect. [3] The distribution of angular momentum values and approxi- 
mations used in the derivation of the angular momentum correlation 
function are discussed in Sects. |4] and [5] respectively. Analytical 
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derivations are contrasted with the results from «-body simulations 
of cosmic structure formation in Sect. |6] and to observations of an- 
gular momentum induced galaxy shape correlations in Sect.|7] Im- 
plications of such spin-induced shape correlations for weak grav- 
itational lensing are derived in Sect. [8] and a summary in Sect. |9] 
concludes the review article. 

The choice of specific values for the cosmological parameters 
is not as relevant for angular momentum correlations compared to 
other cosmological observations, because the angular momentum 
build-up of galaxies takes place at early times when the universe is 
matter dominated. Nevertheless I will attempt to keep the deriva- 
tions as general as possible and use as a framework spatially flat 
dark energy cosmological models with Gaussian adiabatic initial 
perturbations in the cold dark matter density field. 



2 COSMOLOGY AND STRUCTURE FORMATION 

This section provides a compilation of dark energy cosmologies 
(Sect. im , the description of the fluctuations statistics of the den- 
sity field with power spectra (Sect, \2.2i . the structure formation 
equations and their linearisation (Sect. l2.3T ). the linear growth equa- 
tion (Sect. |2.4t . Lagrangian perturbation theory (Sect. |2.5T > and halo 
formation by spherical collapse (Sect. |Z6l ). 



2.1 Dark energy cosmologies 

In a spatially flat dark energy cosmology with the matter density 
parameter Q.,,,, the Hubble function H(a) = dlna/d/ is given by 
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+ (1 - n„)exp 



da' 



1 + w(a') 
a' 



(1) 



with the dark energy equation of state w{a), and omitting the radi- 
ation density fl,-, which has a negligible influence on the late-time 
expansion. The value vv = - 1 corresponds to the cosmological con- 
stant A. The relation between comoving distance x (given in terms 
of the Hubble distance = c/Hg) and scale factor a is given by 



r' , 1 

; da —; , 



(2) 



with the speed of light c. The dark energy equation of state w(a) is 
conveniently approximated by its first order Taylor expansion with 
respect to the scale-factor a, 

w(a) = wo + (l- a)Wa, (3) 

introduced by IChevalUer & Polarskil i200lh and iLinder & JenkinsI 
( 12003 '), for non-interacting dark energy models with a slow time 
evolution w{a). The cosmic time in units of the Hubble time tn = 
1 /Ho follows directly from the definition of the Hubble function, 



: f da' — 



(4) 



For a number of purposes, in particular as the time variable in the 
structure formation equations, the definition of the conformal time 
T] in analogy to the comoving distance x from the differential d;; = 
At la = da/(a^H{a)) is useful. 



2.2 CDM density spectrum 

The fluctuations of the cosmic density field 6(x), which is defined 
as the fractional perturbation of the density field p(x) relative to the 
background density (p) = Q.„,Pcni, 



Six): 



p(x) - (p) 

(P) ' 



(5) 



are assumed to be Gaussian and are described by the variance be- 
tween two Fourier modes 6(k) = j d^x d(x) &xp{-ikx). 



(6{k)6{k')) = {InfSoik + k')P(k), 



(6) 



which defines the power spectrum Pik) in the case of homogeneous 
and isotropic fluctuation statistics. Inflationary models suggest the 
ansatz 

P(k)oc k"'T\k) (7) 

jGuth & Pi|[T982h for the power spectrum with the CDM transfer 
function T( k), which is we l l appr oximated with the polynomial fit 
proposed bv lBardeen et al.l ( Il986h , 



T(q) 



ln(l + 2.34^) 
2.34q 



(l + 3.89^ + (le.lqf + (5A6qf + (6Jlqf) 



The wave vector k = qT is rescaled with the shape parameter F, 
which describes the peak shape of the CDM power spectrum P(k). 
[Sugivama ( 1995) gives a convenient fit to numerical data in cos- 
mologies withflm t 1, 



:n„,.exp(-a,(i.|^]]. 



(8) 



where F is measured in units of (Mpc//;)"', such that q is dimen- 
sionless. The normalisation is given by the variance erg of the den- 
sity field on the scale = 8 Mpc//? at the current cosmic epoch. 



dkk'^P(k)W^(kR), 



(9) 



with a Fourier transformed spherical top hat filter function, W(x) = 
371 (jf)/ X. j((x) is the spherical Bessel function of the first k ind of 
order i jAbramowitz & Stegunlll972l : lArfken & Webej|2005h . 



2.3 Structure formation with CDM 

The formation of cosmic structures is based on gravitational ampli- 
fication of seed perturbations in the cosmic density field. Being a 
hydrodynamical self-gravitating phenomenon, structure formation 
is described in the comoving frame by the system of differential 
equations composed of (;) the continuity equation 



d 

— 5 + div[(l + 



0, 



(10) 



which relates the time-evolution of the density field to the diver- 
gence of the matter fluxes 7 = (1 -1- S)v, (ii) the Euler-equation 



■^v + aHv + uVt! ■ 
drt 



-VO, 



(11) 



which describes the evolution of the peculiar velocity field v from 
the gradient VC) of the peculiar gravitational potential <b, acting on 
a fluid element, and finally {Hi) the Poisson-equation 



AO = -a,„{r]) {aHfS, 



(12) 



which gives the gravitational potential O induced by the matter dis- 
tribution 5 (Newton's constant has been replaced with the definition 
of the critical density, and absorbed into the ambient density into 
Q.m{rj) with the conformal time 77 as the time variable). The expand- 
ing background is described by the Hubble function H(a) and the 
usage of Newtonian dynamics and Newtonian gravity is well jus- 
tified on sub-horizon scales. Due to the fact that cold dark matter 
is coUisionless, viscous and pressure forces are negligible and the 
formati on of objects takes plac e by phase-space mixin g, dynamical 
friction ( IChandrasekhaj[l943h and violent relaxation 1 Lvnden-Belll 
1967). 

Linearisation of the structure formation equations by substi- 
tuting a perturbative expansion of the density- and velocity field 
yields the linearised continuity equation. 



—6 + divt) 
drt 



0, 



(13) 



and the linearised Euler-equation, 

d 

—V + aHv = -VU), (14) 
dl 

which are valid if the overdensity 6 is small, S ^ I. 

It is instructive to decompose the velocity field into the diver- 
gence 6 = diwv and the vorticity = rotv, for which one obtains 
the evolution equations from the linearised Euler-equation 

■^e + aHe+-n^(T])(aHfS = and -^u; + aHw = 0. (15) 
orj 2 oj] 

Solving the latter equation shows that the vorticity decreases 
oc 1/a such that any initial vortical excitation is wiped out 
rapidly, which disfavours models in which galactic angular mo- 
menta ar ise from vortical init i al pe r turbations in the veloc- 
ity field IZ el'dovic h & Novikovl Il970l : IJonesI 1 1 9761: iGotll 1 1 9771; 
lEfstathiou & Silk 1983). Vorticity can then only be generated by 
a nonlinear amplification term of the form iol(v x cj) in the full 
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Euler-equation if there are non-zero stresses. The flows which de- 
velop in a collisionless fluid are necessarily laminar potential flows. 
Only collective processes between the dark matter particles can 
emulate viscosity and pressure, which happens on fully nonlin- 
ear scales, and generate vorticity by the nonlinear te rm rot(u x cj) 
jPichon & Bemardeau|[l99^ : iBemardeau et al.ll2002h . 



2.4 Linear structure growth 

The linearised continuity equation can be combined with the time 
derivative of the evolution equation for the velocity divergence to 
form the growth equation, whose solution describes the homoge- 
neous growth 5{x,a) = D+(a)S{x,a = 1) of the density field in 
the l inear regime of structure formation, \S\ <g: 1 ( |Peebles|[T980l : 
[Turner & White 1997; Under & Jenkins 2003), 



d\nH\ d 



1 , 

—£)+(«) + - 3 + 
da-^ a\ dlna / da 



— = ^£l„(a)0+(fl). (16) 



The growing solution D+(a) of the growth equation's two solu- 
tions is referred to as the growth function and assumes the sim- 
ple solution D+(a) = a in the SCDM cosmology, where Sl„, = 1 
and H(a)/Ho = a"^'-. This scaling motivates the usage of the 
scale factor a as the preferred time variable, and to transform the 
differentials with d/drj = a-H(a)d/da. Hence, the initial condi- 
tions for solving the growth equation are given by D+(0) = and 
dD+(0)/da = 1, due to matter domination in early times. The term 
3 + d In H/d In a is equal to 2 - with the deceleration parameter 
q, and acts as a fiction term which suppresses structure formation 
on scales whose dynamical time scale is larger than the time scale 
of the Hubble expansion. 

Homogeneous structure formation is equivalent to the notion 
of independently growing Fourier modes, 

6(x, a) = D+(a)6(x, a = I) ^ S(k, a) = D+(a)S(k, a = I), (17) 

which conserves the Gaussianity of the initial conditions. The 
Gaussianity of the initial density perturbations is a consequence of 
inflation, where a large number of uncorrelated quantum flucuta- 
tions are superimposed, yielding a Gaussian amplitude distribution 
due to the central limit theorem. As long as the growth of struc- 
ture is linear and thus homogeneous as described by eqn. (|T7]l, all 
modes growth independently conserving the Gaussianity of the ini- 
tial density field. 



to the displacement potential 4' and allows to write by using the 
Poisson equation 

' d2 



J(q, ri) div 



^ ^ d 
dr]^ dr] 



(20) 



with the mass conservation condition (1 + (5)d^x = d^q. Therefore, 
1 -I- S(x, rj) - l/J(q, rf) with the Jacobian 

J{q, n)='^ = (det [S,j + didj^]y' , (21) 

meaning that the tidal forces didjV cause a compression of the cos- 
mic density field. The limit of applicability of Lagrangian pertur- 
bation theory is reached if 7 = 0, because the density field for- 
mally diverges. This happens if two trajectories cross and if the 
mapping q ^ x becomes ambiguous. The linear solution to 
can be derived as t^iq, rj) = D+(r])S{q), using the solution D+(r]) 
to the homogeneous growth equation (Catelan 1995). With this 
potential and the perturbative mapping x(q, q) = q — V^{q,ri) 
is refered to as the Zel'dovich-approximation jZel'dovichl 1 197d : 
lBucherl|[l992l) . The analytical power in the Zel'dovich approxima- 
tion lies in the extrapolation of nonlinear dynamics from the initial 
conditions, without having to follow the evo lution of nonlineari- 
ties computationally in t he Eulerian approach jBouchet et al.ll992 ; 

ISahni&Colesil995 : 



ties computationally m t he buienan approach 
Buchert & Ehlersll 19931: Isouchet et aij|l995l :[ 



Catelan et alJll995l ; lBemardeau et alj2o'm 



2.6 Galaxy formation and biasing 

Galaxies form at local peaks in the large scale structure by ellip- 
soidal collapse (Sheth et al. 2001) of the dark matter perturbation. 
In the collapse, the halo decouples from the local Hubble expan- 
sion which is overcome by the halos own gravitational field. The 
baryonic component inside a dark matter halo undergoes radiative 
cooling and forms a rotationally supported galactic disk. 

The fractional perturbation A«/;2 in the mean number density 
n of dark matter haloes hosting a galaxy is related to the overdensity 
field S = Ap/(p) by the bias parameter b, 



An 



(22) 



which is in general scale dependent I Lumsden et al. Jl989l) and 
slowly time-evolving ( lFrvlll996l:lTegmark & Peebleslll998h . 



2.5 Lagrangian perturbation tlieory and nonlinear growtli 



3 TIDAL TORQUING 



The appropriate tool for describing the angular momentum build- 
up of haloes in the large-scale structure, which is a perturbative 
pro cess, is Lag rangian perturbation theory (for a n detailed review, 
see lSahni &"Co les 1995; Bernarde au et alj|2002h . The central ob- 
jects in Lagangian perturbation theory are the particle trajectories 
linking the initial positions q to their positions x at time q, 



x(q,q) = q-W{q,q), 



(18) 



with the displacement potential T which describes the potential 
flows developing in the large-scale structure. The equation of mo- 
tion for a particle is given by 

d2 



-X + aH — X 

dq'- dq 



-VO. 



(19) 



The divergence of this equation yields the connection to the Eu- 
lerian approach by relating the physical gravitational potential <b 



This section treats tidal torquing as a mechanism for angular mo- 
mentum build-up of dark matter haloes in detail. After a short 
qualitative description (Sect. lSTTT l. I introduce tidal torquing in La- 
grangian perturbation theory using the Zel'dovich approximation 
as a dynamical model (Sect. [T2] l and discuss the time evolution of 
the angular momentum for constant torques in various cosmologies 
(Sect. 13. 3t . A crucial technical point is the misalignment between 
the eigensystems of the inertia tensor and the tidal shear, which I 
elaborate on in Sect. l3.4l followed by a derivation of the scaling of 
the angular momentum with halo mass (Sect. [331 l. 

3.1 Mechanism of tidal torquing 

The current paradigm for generating the angular momentum of 
galaxies is tidal torquing, where the tidal gravitational fields exerts 
a torquing moment on the protogalactic object prior to collapse. It 
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should be emphasised that in this picture the angular momentum 
does not originate from the vorticity which is amplified by the 
nonlinear term iot(v x u>) in the Euler-equation as explained in 
Sect. 12.31 but is rather generated from (vorticity-free) shear flows, 
in which protogalactic objects are embedded. During the shearing, 
each protohalo is deformed and acquires a rotational motion com- 
ponent. Eventually it decouples from the shear flow by collaps- 
ing under its own gravity, which reduces the length of the lever 
arms and makes torquing inefficient. For the perturbative descrip- 
tion of the deformation of the object by the action of shear flows, 
Lagrangian perturbation theory is suitable. The difficulty in this ap- 
proach lies in the fact that the shear flows themselves influence the 
boundary surface of the embedded object and hence its total mass 
and inertia. Tidal torquing effectively solves the problem of gen- 
erating vorticity from a laminar flow by introducing a non-simply 
connected density and velocity fields, because during spherical col- 
lapse the protohalo volume is separated from the ambient fields. 

This me c hanism was first q uant itatively invest igated by 
iDoroshkevicl] h970h . lwhit3 ill984) and lWessod ( Il985l) . building 
on the original idea by HovlellT949l) and Sciamal jl955l) . Assum- 
ing a non-spherical shape of the protogalactic region, the angular 
momentum grows at first order and linearly in time in Einstein- 
de Sitter universes, whereas in spherical regions, the acquisition of 
angular momentum is only a second order effect due to convec- 
tive m atter streams on the boundary surface, as shown bv |Peeble3 
( Il969l) . 



3.2 Acquisition of angular momentum by tidal torquing 

Quite generally, the angular momentum L of a rotating mass distri- 
bution p(r, t) contained in the physical volume V is given by: 



Euler frame 



Lagrange frame 



m 



I'' 

Jv 



r (r- r)x v(r, t)p{r, t). 



(23) 



where v(r,t) is the (rotational) velocity of the fluid element with 
density p(r, t) = (p)(l + 6(r, t)) at position r around the centre of 
gravity f. In perturbation theory, 6 <s: I, and the density field can be 
approximated by assumi ng a constant density (p) = flmPcrit insid e 
the protogalactic region. JWhite ( 1984'),'Catelan c& Theun^ l ll996ah . 
iTheuns & CatelM] ll99'A) and Crittenden et al. (2001) have de- 
scribed the growth of perturbations on an expanding background 
in Lagrangian perturbation theory: The trajectory of dark mat- 
ter particles in comoving coordinates to first order i s given by 
the Z el'dovich approximation as the dynamical model l lZerdovichI 
Il970 ): 

x{q, t) = q- D+(t)W'i'(q) x = -D+VP, (24) 

which relates the initial particle positions q to the positions x at 
time t. The particle velocity x follows from the Zel'dovich-relation 
by differentiation by the time- variable. The growth function Z)+(f) 
describes the homogeneous time evolution of the displacement field 
and contains the influence of the particular dark energy model. 
In the Lagrangian frame, the expression for the angular momentum 
becomes 

L = poa^ I d^^ (x - x) XX =^ poci^ I d^q (q - q) X x, (25) 

Jvi Jvi 

where the integration volume is defined in comoving coordinates 
as well. Assuming that the gradient Wiq) of the displacement 
field ^'{q) does not vary much across the Lagrangian volume Vt, 
a second-order Taylor expansion in the vicinity of the centre of 
gravity q is applicable: 





Figure 1. Principle of tidal torquing: A linear variation of the displacements 
across the protogalactic object in the space-fixed Euler-frame translates to 
a rotational motion in the comoving Lagrange-frame after collapse. 



d^^iq) - 3„>P(9) + Yj(q - 9)/iT„/,, 



(26) 



The expansion coefficient is the tidal shear at the point q: 

'f'.riq) = d^dy^iq), (27) 

because the Zel'dovich displacement field T is related to gravita- 
tional potential O and can be computed as the solution to Poisson's 
equation AT = 6 from the cosmological density field d. The gradi- 
ent d„^(q) of the Zel'dovich potential displaces the protogalactic 
object, which is neglected in the further derivation, as I only trace 
differential advection velocities responsible for inducing rotation. 
Identifying the tensor of second moments of the mass distribution 
of the protogalactic object as the inertia 1^^, 



I d^'^ (9 - q)i3(q - q)a 



one obtains the final expression of the angular momentum L„ 



(28) 



(29) 



Here, eqn. l |29l l implicitly assumes a peak constraint because the 
inertia of the protogalactic volume is needed, which can only be 
defined in a sensible way for a density peak from which a halo 
forms by collapse. 

Fig. [T] illustrates the torquing action on a protogalactic region 
in the Euler- and Lagrange frames: If the Zel'dovich displacements 

vary linearly across the protogalactic cloud, they introduce a 
rotational motion in the comoving Lagrangian frame. Being pro- 
portional to the second derivative of the potential T, the rotation 
is determined by the tidal forces "P^^ = (9^5y*P. Tidal torquing is 
effective until the moment of turn-around in the spherical collapse 
picture, because the collapse dramatically reduces the lever arms. 
After the collapse, the halo conserves the angular momentum it has 
accu mulated until turn-arou nd . 

IPorciani et al.l ( l2002al) give an intuitive argument on the ori- 
entation of the angular momentum relative to the principal axis 
system of the tidal shear, which is a very useful relation: In the 
eigenframe of the tidal shear, the components of the angular mo- 
mentum are simply given by L„ oc IpyCVp — "Py), where the indices 
are cyclic permutations of (1,2, 3). "P„ with a single index are the 
three eigenvalues of T. Averaging over all orthogonal transforma- 
tions relating the eigensystems of I and 1* yields the result that the 
largest component of L would be the one for which I^P^ - "P^l is 
largest. Because of the ordering < ^2 < 4^3, this is necessar- 
ily L2, Li oc I^Pj - Til, i.e. the angular momentum tends to align 
itself parallel to *P2, i.e. the axes corresponding to the intermediate 
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Figure 2. The influence of the dark energy model on the time evolution of 
angular moinenta: Q(a) as a function of scale factor a for ACDM (wq = - 1 
and Wa = 0, solid line), for quintessence (wq = -2/3 and Wa = 0, dashed 
line) and a dark energy model with variable equation of state (wq = -2/3 
and Wa = 1/3, dash-dotted line). The matter density parameter was chosen 
to be f2„, = 0.25 in all cases. 



eigenvalue of the tidal shear. This argument is not generalisable to 
the case with correlations between shear and inertia: In the ideal 
case of a Gaussian random field, one can transform into the eigen- 
system of the shear 1*, such that only the off-diagonal components 
of I are used in the computation of L, which are statistically un- 
correlated, and this argument fails if there are contributions from 
the diagonal components of I, which can disturb the ordering of the 
angular momentum components. 



3.3 Time evolution of the angular momentum 

For investigating the influence of dark energy on the spin-up via the 
Z)+ -factor, it is convenient to rewrite the time dependence of Z)+ in 
terms of the scale factor a by dZ)+/d? = aH(a)dD+/da, yielding: 



3 dD+ 

La = a H{a) — 6„^y ^ Ipa-^a 



(30) 



Fig. |2] compares the time evolution of the angular momentum in 
dark energy cosmologies with SCDM. I define the ratio 



with 



gDE(Q) 
<?SCDM(a) 



q(a) = a^H(a) ^^* , 
da 



(31) 



(32) 



where D+{a) is normalised to unity today. In SCDM these for- 



mulae simplify to H(a)/Ho 



,,-3/2 



D+{a) = a and consequently 



?scDM = Hqu^^-^. Fig.|2]compares Q(a) for dark energy models with 
a constant and a time-evolving equation of state to a model with 
a cosmological constant A and suggests that the spin-up of haloes 
in dark energy models half as fast compared to SCDM, and the 
choice of the equatio n of state affects the time evolution signifi- 
cantly ( ISchafeJl2008l) . 



3.4 Misalignment of the shear and inertia eigensystems 

In order to illustrate the consequence of the contraction of the prod- 
uct X = IT with the Levi-Civita-symbol e„py in the definition of the 
angular momentum, one can carry out a decomposition of this ten- 
sor, X = X^+X" (Schafer 2008), into an antisymmetric contribution 
X", defined via the commutator [I, T], 



X, 



(33) 



with the symmetry (X")' = ^(IT - Tl)' = |(TI - \W) = -X under 
matrix transposition (I and T are symmetric matrices) and into the 
corresponding symmetric matrix X^ by using the anticommutator 
jl, *Pj between inertia I and tidal shear 'V: 



(34) 



with (X^)' = +X^. In the derivation of the angular momentum L, 

(35) 

cr 

the permutation symbol e^py picks out the antisymmetric contribu- 
tion X", by virtue of e„py{X^^ + X^^) = X^^, because the contrac- 
tion of the symmetric tensor X* with the antisymmetric permuta- 
tion symbol vanishes, e^pyX^^ = 0. 

Hence, the protogalactic objects will only acquire angular mo- 
mentum if the commutator X" between the inertia and the tidal 
shear is non-zero, which means that the inertia and shear tensors 
are not supposed to be simultaneously diagonisable, i.e. they are 
not allowed to have a common eigensystem and have to be skewed 
relative to each other. An extreme example of the dependence of 
the angular momentum on the commutator X" is a spherical pertur- 
bation, which will never start rotating irrespective of the strength of 
the tidal fields because one can always choose the eigensystem of 
the inertia to coincide with that of the tidal shear, resulti ng in a van- 
ishing commutator. In order to capture this mechanism. iLee & PenI 
|2000) and Crittenden et al. (2001) have used an efi"ective, param- 
eterised description of the average misalignment of the shear and 
inertia eigensystems, gauged with numerical n-body data. 

The symmetric contribution X^, which measures the degree 
of alignment of the inertia and shear eigensystems, causes an 
anisotropic deformation of the protogalactic region during the 
course of linear structure formation prior to gravitational collapse. 
Consequently, the determination of ellipticity distributions is likely 
to be affected even in the stage of linear structure formation, and 
predic tions of triaxiality ba sed on peak shapes in Gaussian random 
fields teardeen et al.|[l98^ could be refined using an adaptation 
of the formalism outlined above. Additionally, the tidal forces de- 
scribed by X* have an influence on the total mass of the object 
whi ch will form by ellipsoidal collapse on a particular peak (see 
e.g. iLee et ai]|2005h . In a very detailed studv. lDubinskil ( 1 19921) in- 
vestigated the impact on tidal shearing on halo formation in an n- 
body simulation. He found significant deviations from the shapes 
of peaks in the initial density field compared to the halo triaxiality 
distribution, by the action of X^, and influence of the tidal forces 
on the boundary defining the dark matter object which forms at a 
peak in the density field. 



3.5 Scaling of angular momentum with halo mass 

In the CDM paradigm, dark matter haloes form by spherical col- 
lapse, which decouples the dark matter object from the tidal forces. 
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Hence, the angular momentum built up until the moment of turn 
around (at which the radius stalls) can be expected to be a good es- 
timate of the total anaular momentum. As shown bvl Peeblesl i fT96^ 
for haloes in an SCDM cosmology, more massive objects have 
larger angular momenta, the scaling being given by L oc M^'^ . In a 
low-density cosmology like ACDM, the overdensity at turn around 
is only slightly larger than unity, 6 = 1.07, such that the continu- 
ity equation can be written as 5 = -D+A'F, which is equivale nt to 
setting P+ ^ 1/Ay at turn around. ICatelan & Theunj jl996bh and 
ISugerman et al.l ( I2OO0I) continue by denoting the total mass of the 
protohalo by M and its radius at turn around by R, and deduct that 
at turn around 



L - a'D.A^ MR- - £lp-y^M''' cc O-""^ 



(36) 



which gives the scaling L oc M^'^ for the relation between angular 
momentum and halo mass. The term MR^ is an estimate of the 
inertia of the protohalo. In summary, the origin of the mass-angular 
momentum relation is the fact that the angular momentum can grow 
until the time of turn-around, which is determined by the halo mass. 



4 ANGULAR MOMENTUM DISTRIBUTION 

This section deals with the statistics of the angular momentum dis- 
tribution. I provide a formal derivation of the moments of the an- 
gular momentum distribution from a Gaussian random process, re- 
stricted to galaxy formation sites (Sect. l4Tt . followed by a deriva- 
tion of the angular momentum variance, which uses the distribu- 
tion of protohalo shapes in a Gaussian random field as an input 
for the inertia distribution (Sect. l42t . The restriction of galaxy for- 
mation to local maxima in the density field results in biased, non- 
Gaussian angular momentum distributions, which are discussed in 
detail (Sect. I4.3t . Finally, the angular momentum statistics is ex- 
tended to nonlinear structure formation (Sect. I4.4K followed by a 
compilation of properties of the spin-parameter A and its influence 
on disk formation (Sect. l43t . 



4.1 Moments of the angular momentum distribution 

The probability density p(L)dL of the galactic angular momenta L 
follows from a joint Gaussian random process, which returns the 
variables which are necessary to describe the angular momentum 
acquisition, restricted to local peaks in the large-scale structure be- 
cause galaxy formation is thought to take place at local peaks only. 
These variables are density field S(x), the density gradient da{x), 
the second derivatives 6„n(x) and the tidal field ^^aix): 



r d^fc 
J W 



^(^) = I 7:r^<5(*)exp(iAjc), 

Sa(x) = — — =1 ——rk„S{k) expiikx), 

dx„ J (2ny 



d^6(x) 



J (2/r)3'" 



kpSik) exp(iAjt:), 



(37) 
(38) 
(39) 



The tidal shear follows from the solution of the Poisson equation 
Al'(x) = 6{x) linking the Zel'dovich displacement potential ^{x) 
to the density field S{x): 



"V^pix) 



d^^{x) r d^k k„k„ 



dxadxp 



(40) 



An important consequence of eqns. l |39t and ( I40t will be the fact 
that the angular momentum correlation is determined by two mech- 
anisms with differing correlation length: a short range correlation 
of the peak shapes and hence the inertia, and a long range corre- 
lation mediated by the tidal shear. There is an interesting analogy 
to solid state physics: Electron spins in ferrimagnetic materials are 
subjected to interactions of two different correlations lengths, too. 
The occurrence of two correlation lengths is crucial to the angular 
momentum generation. The longer correlation length of the tidal 
shear ensures that the tidal field has contributions from the ambient 
density field, and is not sourced by the perturbation alone. If the 
latter were the case, shear and inertia would be degenerate and the 
protohalo could not acquire angular momentum. 

A peculiarity worthwhile mentioning is the fact that the 
density field S(x) is degenerate with the trace ir^apix) of the 
tidal shear beca use of Poisson's equation AT = trTj,^ = 6 
jHeavens & Pea cock 1988). For that reason, the density field will 
appear in the above outlined random process as a derived quantity, 
whereas the entries of the tidal shear matrix will be drawn from the 
Gaussian distribution, with the peak restriction in place. As a conse- 
quence, the multivariate Gaussian distribution has only 15 degrees 
of freedom, instead of 16. A second reason why this degeneracy is 
sensible is the fact that the trace of the tidal shear, being a coordi- 
nate system independent quantity, transforms as a scalar just like 
the density field itself. 

The peak restriction is modeled by the criterion that the Gaus- 
sian process is required to exceed a threshold v in density which 
corresponds to a minimal peak height for galaxy formation, 6{x) > 
V, to be of vanishing density gradient, 6a{x) = 0, and of negative 
curvature, Sap{x) < 0. In addition, the curvature tensor 5„^(x) < 
is used to derive the inertia of the protogalactic region. 

The number density of maxima in the density field can be 
modeled by a point-process 1 Bardeen et al.lll985 iRegos & Szalavl 
ll995l : lHeavens & Shethlll99sl) : 



Y^^lix-x-^ 



Close to the maximum i at Xi 
gradient 6a(x) is applicable: 

Saix) = 2 5ap{Xi){X - Xi)ii. 



(41) 



a Taylor expansion of the density 



(42) 



With this expansion on obtains for the peak density: 
"peak = '^6%{Slp(Xi)8a') = |det(5„^| i5^(<5„). 



(43) 



The relation takes a simpler shape when considering the eigensys- 
tem of the mass tensor -Sap'- Being a symmetric tensor, it has the 
three real eigenvalues /i,, i = I ... 3, which allows to replace the 
determinant by the invariant quantity |/li/l2/i3|. 

The constraints can be combined in a mask C{v), which is de- 
fined as a function on the vector v containing the derivates of the 
Gaussian random field under consideration: Peaks in the density 
field are defined as points with amplitudes in units of the variance 
o"o = (<5^) exceeding a certain threshold v and exhibiting a vanish- 
ing gradient Sa as well as negative curvature Sap'- 



C(u) = Si [Saix)] \AiA2Ai\ @{Ai) [Six) - crov] . 



(44) 



0(x) denotes the Heaviside step-function. By arranging all random 
derivates in a vector v, the peak density Wpeak, i-e. the expecation 
value for the number density of peaks in the fluctuating density 
field 5 which exceed a threshold viTo can then be derived from the 
multivariate Gaussian random process p{v)dv. 
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^peak 



dv p{v)C{v), 



(45) <L2> = ^a'bl - {I„„f) j ^ k^'P^ 



which corresponds t o the integral of the differential peak density 
"peak(i')dv defined bv lBardeen et alj ( Il986l) . 

In analogy, the moments of the angular momentum distribu- 
tion with a restriction to peak regions in the density field can be 
obtained with (Schafer 2008): 



<L") = \— { dv p(v)C(v)L''(v). 



(46) 



Under the assumption of isotropy, it is sufficient to consider the 
distribution p(L)dL of the angular momentum along a single di- 
mension. In both formulce, the normalisation factor 1/npcak ac- 
counts for the discreteness of the measured quantity. In principle, 
all peak restricted moments (L") can be derived and the distri- 
bution p(L )AL reconstructed using the characteristic function ip{t) 
jArfken & Weber 2005i) . 



j dLp(L)exp{-iLt) = J]{L"}'^, 



(47) 



with <L") = / dL L"p(L), and successive inverse Fourier trans- 
form. The function C(v) makes the random process effectively non- 
Gaussian and breaks the property of Gaussian distributions that the 
even moments are powers of the variance cr^. 

A very important result is the anticorrelation of the angular 
momentum magnitude with the peak height v, i.e. especially the an- 
gular momentum variance decreases if peaks of higher overdensity 
are selected. iBj umenthal et al. ( 1984), Hoffman ( 1986), Hoffman 
( Il988l) and lHeavens & Peacock! ( Il988h investigate this effect in an 
analytical derivation as well as by computing a realisation of a 
Gaussian random field and measuring the tidal shear and inertia 
at peaks, deriving results on the statistics of the angular momen- 
tum and of the spin parameter /i as a function of v = S/ctq: The 
spin of an object depends on the time during which its dynamics 
is described by linear perturbation theory and during which tidal 
torquing is effective, and on the correlation of inertia and tidal shear 
tensors, which depends on the randomness of the field. Typical vari- 
ations of the angular momentum amount to a decrease of ~ 30% in 
amplitude if the peak height is increased from y = 1 to v = 2. 



4.2 Angular momenta and peak shape distributions 

ICatelan & TheunsI ( Il996al) carry out a computation of the variance 
(L^) resulting from an unconstrained random process, effectively 
by embedding ellipsoidal haloes with axis ratios drawn from the 
correct multivariate Gaussian distribution at random positions in 
the shear flows of the large-scale structure, and generalise this to the 
derivation of the angular momentum probability distributions. They 
assume no correlation between the tidal shear and the inertia and 
start with an evaluation of the tidal shear variance (using statistical 
isotropy) 



(48) 



(49) 



Pa-^ P' a-') = (Spa-Sp'a-' + Spp'Saa' + Sp^r'Spia^ O, 

with the abbrevitation 

^ji r dk 
I* = T7 7T^ k''P^(k)W\kR), 
15 J (Iny 

directly computed from the spectrum P'^ik) of the Zel'dovich- 
potential "P, smoothed on the scale R. In this way, they can link 
the angular momentum variance to the shape of the protogalactic 
object, described by its inertia 



{k)W\kR). (50) 



This expression, being scalar and independent of the coordinate 
system, allows to transform into the principal axis system of the 
inertia tensor and to simplfy the expression. 



• (/„„)- = 2{pi\ - 3p2), 



(51) 



where pi and p2 are the invariants of the inertia tensor: pi = tri = 
Li + i2 + i3, and p2 = tit2 + L1L3 + 1213, formed from the eigenvalues 
i, , i = 1, 2, 3, of the inertia tensor. Specifically, the quantities pi and 
P2 are invariant under orthogonal transformations, and in addition 
there exists a third invariant, which is the determinant p^ = 111213 
of the inertia tensor. A spherical object has 3 identical eigenvalues 
ii = I2 = 13 = i, and therefore pi = 3l, p2 = 3r, and hence 
the expression p^^ - 3p2 vanishes, in accordance with the fact that 
spherical protohaloes can not acquire angular momentum by tidal 
torquing. For a non-spherical volume, p^^ - 3p2 > 0, such that tidal 

torq uing can be effectiv e. 

ICatelan & TheunsI l ll996ah now link the angular momentum 
distribution to the distribution of peak shapes by estimating the in- 
ertia of a halo forming at a peak from the local curvature of the 
density field: Clearly, a small value of the local curvature corre- 
sponds to a large object, with a large inertia. They show that the 
term yu^ - 3p2 replacing the inertia in eqn. l lSOt depends on the el- 
lipticity e and prolateness p of the protogalactic object. 



p] - 3p2 



2"3'' 2 2/^o"oV Me,p) 



52 



(52) 



^xo-2/ SHe,py 

with 7/0 = p,„a'. The ellipticity e and the prolateness p are related 
to the eigenvalues /!, , ; = 1, 2, 3, of the curvature tensor dad/36. 



P 



2(^1-1-^2 + ^3)' 

A[ — 2.A2 + /I3 



(53) 



(54) 



2(Ai +A2 + Ai) 

For a triaxial ellipsoid, the parameter e > measures the ellipticity, 
and p gives information whether the object is oblate (0 < p < e) 
or prolate (-e < p < 0). For an alternative description of the local 
curvature in terms of e and p instead of the eigenvalues Aj, i = 
1, 2, 3, the additional parameter x, 

A, + Al + A3 

X = ^ (55) 

0-9 

is needed. The moments cr„(R) derived from the power spectrum of 
the smoothed potential can be computed using 



P>i(k)W\kR). 



9 r dk 

The polynomials ^ and S are given by: 

:n(e, p) = [p{\+ p)f + 3e^{l -6p + 2p^ + 3e^), 
S(e,p) = (l-2p)[(l+pf-9e^]. 

Using these relations, the final step is to relate the angular mo- 
mentum distribution to the proba bility distribut i on of peak shapes 
p{v, X, e, p)dv dx de dp derived in lBardeen et al.l jl986t) . 



(56) 

(57) 
(58) 



(59) 



p(y, X, e, p) oc x^-Wie, p) exp (~ - ^-x\3e' + p'-) - %t^. 

\ 2 2 2(\ -y'- 

with the function "Wie, p) = e(e^ - p^)S(e, p). The probability dis- 
tribution is defined in the allowed triangular domain < e < 1/4, 
-e < p < e and 1/4 < e < 1/2, 3e - 1 < p < e. In a quantita- 
tive investigation. ICatelan & TheunsI ( fl996al) find that by choosing 
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typical peak heights, the value of the specific angular momentum 
of spiral galaxies can be reproduced, but they find discrepancies 
for higher peaks corresponding to elliptical galaxies, as the specific 
angular momentum predicted is too large. Another important result 
that the median of derived distribution of the spin parameter A from 
this model agrees well with the value observed in spiral galaxies. 



ing again the Poisson equation, which relates the density field to 
the trace of the tidal shear. 

In particular, this fortunate coincidence can not be used in the 
case of angular momentum correlations, for which the full correla- 
tions between the inertia and shear across two points in the large- 
scale structure as well as between themselves at the same location 
has to be taken account of, as will be demonstrated in Sect. 15. II 



4.3 Angular momentum biasing 

The restriction of the random process for angular momenta to 
peak locations with has two consequences, which effectively in- 
troduce a biased angular momentum distribution. Firstly, as shown 
by [Heavens & Peacock ( 1988), who perform a direct computation 
of the distribution p{L)AL by setting up the joint Gaussian ran- 
dom process outlined above and integrating out the peak restric- 
tion in the derivation of p(L)dL, the angular momentum distribu- 
tion is non-Gaussian and although the mean angular momentum 
needs necessarily to be zero due to isotropy (a nonzero mean an- 
gular momentum of the large-scale structure would define a prefer- 
ential direction), the odd higher order moments of the angular mo- 
mentum distribution have no reason to vanish. Additionally, they 
provide fitting formulse for the resulting angular momentum distri- 
bution p{L)AL. 

The second consequence of the peak restriction is a reduced 
variance of the angular momentum field, compared to that of a con- 
tinuous field, due to the covariance of the curvature of the density 
fie ld with the tidal shear. The magnitude of this effect is determined 
bv lCatelan & Theunj ( [T9963) . who show that the restriction to peak 
location reduces the variance of the tidal shear by (1 - y^). 



with the value y, 
7 



1 



(60) 



(61) 



This reduction, which is mass-scale dependent due to the smooth- 
ing scale R, is due to the fact that the random field is constrained to 
have a peak with a specified shape, and is restricted in its fluctua- 
tions. cr„ are given by eqn. J56b . 

The idea of the derivation outlined in Sect. 14.21 was to place 
protohaloes with axis ratios from the multivariate distribution 
p(v, X, e, p) at random positions in the large-scale structure and to 
torque them up with the local tidal fields, such that one obtains 
the probability distribution of the angular momentum. This proce- 
dure assumes the validity of two simplifications: The fluctuation 
amplitude of the tidal shear field is the same restricted to peak po- 
sitions and in the case of a continuous field, as discussed above. 
Secondly, there are no correlations between the inertia of the proto- 
haloes which form at peak positions and the tidal shear field. This 
is seemingly inconsistent, as the shape of the potential and hence 
the density field needs to be compatible with the shape of the object 
which is forming at a particular position - and one would need to 
take account of the correlations between the inertia and tidal shear 
tensors. In the case of a Gaussian random field, however, it is pos- 
sible to transform in the eigenframe of the inertia tensor, such that 
only the off-diagonal elements of the tidal shear tensor are needed 
for the computation of the angular momentum, due to the occur- 
rence of the Levi-Civita-symbol in eqn. M9\ . These elements turn 
out to be statistically uncorrelated with the diagonal entries of the 
inertia tensor, such that the assumption of uncorrelated tidal shear 
and inertia distributions is justified, the physical reason of this be- 



4.4 Nonlinear tidal torquing 

The theory of angular momentum acquisition by tidal shearing 
has been extende d to nonlinear stag es by using second order per- 
turbation theory jCatelan & Theuns| [l996b) and to include efl^ect s 
of non-Gaussian i nitial p erturbations CCatelan & Theun j Il997h . 
ICatelan & Theun j ( Il996bl) present a perturbative analysis which 
uses higher-order correction terms in Lagrangian perturbation the- 
ory. They find that the variance (L^> of the angular momentum is 
increases by a factor of 1.6 compared to the first-order computa- 
tion, or equivalently, that the angular momenta are higher by about 
30%, and emphasise that the initial torque is a good estimate for the 
total torque acting on the protogalactic object. The second paper 
(Catel an & Theuns 1997) investigates the sensitivity of the angular 
momentum statistics on non-Gaussian initial conditions, using the 
-model and a model with a log-normal distribution of the den- 
sity amplitudes and following tidal torquing in second-order La- 
grangian perturbation theory. Their derivation is generic such that 
it can be applied to both primordial non-Gaussianities in the density 
field or to Gaussianities generated by nonlinear structure formation. 
They find that angular momenta grow from non-Gaussian pertur- 
bations at a lower rate (oc fi^^) compared to Gaussian perturbations 
in an SCDM-cosmology with Q.,„ = 1. The absolute mag- 
nitude of the correctional terms at turn-around can be larger that the 
linear estimates, however, depending on the total skewness. 



4.5 Spin parameter of galaxies and galaxy formation 

Not too surprisingly, the halo angular momentum has far-reaching 
consequences for the internal structure of the galaxy which is 
hosted by a dark matter structure. This is commonly expressed by 
the dimensionless spin parameter A. introduced bv iPeeblesi t l969.) , 

'-gW^^ '''' 

for a halo with angular momentum L, gravitational binding energy 
E and mass M. G denotes Newton's gravitational constant. A cor- 
responds to the ratio between the observed angular velocity of a 
galaxy oj and the angular velocity needed for rotational support ojq: 

(A) L/(MR^) 



A = — 



(63) 



and numerical simulations suggest that the scale radius of expo- 
nential disks is in fact cc A. A depends on the mass scale of the dark 
matter halo, assuming values of up to /I = 0.5 for spirals, and val- 
ues of /t - 0.05 for elliptical galaxies, with a very larg e dispersion 
tefstathiou & Jone j[l979l : Ib ames &Efstathioulll987l) . suggesting 
that dark matter haloes are stabilised mainly by anisotropic veloc- 
ity dispersion. 

The details of the angular mome ntum profile inside a dark 
matter halo is investigated in detail bv lBullocke7ai] ( l200lh . who 
report that the cumulative angular momentum profile M(< j) of 
haloes forming in an n-body simulation is fitted by a function of 
the type 
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M(<7)oc-^, (64) 
J + Jo 

where j = L/M is the specific angular momentum. The profile has 
an implicit maximum value of = jo/(M ~ 1)' is of approxi- 
mate power-law shape for j i jo anf flattens out for j i jo- The 
two parameters yu and 70 are adapted to fit experimental data. The 
angular momentum profile clearly deviate from solid-body rotation, 
implying that dark matter haloes exhibit angular momentum struc- 
ture. The averaged angular momentum directions in shells, how- 
ever, were found to be very w ell aligned with each other. 

Analytical derivations (Peebles '1969'; Heavens & PeacocK 

1 1988.: Steinmetz & Bartelmann 1995; Catelan & Theuns 1996a)^s 
well a s numerical simulations dWarren et alJl 19921 ; Ivitvitska et al.l 
I2OO2I ; iD'Onghia&Navarrol l2007h using realistic CDM transfer 
functions for the initial density field find A to be approximately 
log-normal distributed, 

p(A)dA = -—- exp - (65) 

with the parameters 0.03 ^ p,i < 0.05 and 0.5 < cr,, < 0.7. The 
mean of the distribution is related to /jj by (A) ^ 1.078/J,i for tr^ = 
0.6. 

Using the results from fi-body simulations, IMo et al] j 19981) 
and lAvila-Reese et all ( Il998l) consider the formation of galactic 
disks inside dark matter haloes, while assuming that the baryonic 
component is torqued up to the same specific angular momentum 
as the dark matter, and forms a stable, rotationally supported disk 
with exponential surface density profile, whose mass makes up a 
few percent of the host halo. In this way they are able to repro- 
duce the relation between disk size and rotation velocity, as well as 
the slope and scatter of the Tully-Fisher relation (i.e. the relation 
between rotational velocity and luminosity), calibrating the zero- 
point of the Tully-Fisher relation by fixing the mass to light ratio. 
Further results include rotation curves for galactic disks embedded 
into haloes for a given spin parameter A, which appear very con- 
vincing. The interplay betwe en dark matte r and b aryo nic dynamics 
is ana lysed in the papers by i Tonini et alj j2006l) and iTonini et al.l 
( l2006l) . who find higher spin parameters including baryonic physics 
compared to the dark matter case and a flattening of the central halo 
density profile. 

Other intere s ting results include the paper by 
lEisenstein & Loebl jl995l) . who suggest that quasar progeni- 
tors form by gravitational collapse of perturbations with low 
specific angular momentum on the mass scale of IO^Mq/A. Inside 
these dark matter haloes, they consider the formation of a galactic 
disk on the viscous time scale and conclude that the number 
density of these objects is - 10^^ /(Mpc/h)^ , in accordance with 
the s patial density of quasars. 

Ivan den Bosch et d]( l2002h investigate the spin-parameter dis- 
tribution of the baryonic and the dark matter components and find 
that they are very similar, suggesting that angular momentum is 
conserved during disk formation. While the distributions are simi- 
lar, the angular momentum directions were found to be different by 
a median misalignment angle as large as 30° . In a careful analysis 
of numerical data they conclude that in the course of galaxy forma- 
tion (0 the specific angular momentum of the galactic disk results 
from tidal torquing, (ii) the specific angular momentum of the bary- 
onic component is equal to that of the dark matter, but that (Hi) the 
specific angular momentum is not conserved during disk formation, 
due to virialisation processes in which the baryonic angular mo- 
mentum gets redistributed, while different relaxation mechanisms 
are accessible to the baryonic component and to the dark matter. 



iBailin et alj ( l2005h analyse a set of hydrodynamical simula- 
tions of galaxy formation and find that the disk has a minor influ- 
ence on the overall halo structure, and that the disk is aligned with 
the host halo inside 0. Iryi, , which is a reassuring result for the usage 
of angular-momentum based elli pticity correlation m odels in grav- 
itational lensing. The review bv iMaver et al. I ( l2008h on computer 
simulations of galaxy formation lists some of the numerical issues 
related to angular momenta, e.g. the fact that angular momentum is 
lost artificially by the simulation codes. In addition, they report that 
the angular momentum of a simulated galaxy is resolution depen- 
dent: Increasing the number of particles by a factor of 400 yields 
angular momenta larger by 30%, emphasising the difficulty of nu- 
merical simulations of galaxy formation, especially regarding the 
physics of the interstellar medium. The reliability of the codes for 
determining the alignment of the baryons inside the dark matter 
host structure is still an open question. 

The measurement of the spi n parameter from obser vational 
data is the topic of the papers bv [Hernandez et aL 112003) and by 
lBertaetalJ ( l2008 1. the major issue being that none of the quantities 
in the definition of A (c.f. eqn. I62t is an observable. Both papers 
use SDSS-data and infer A by constructing a model of the galaxy 
consisting of an exponential disk placed in a truncated isother- 
mal sphere, which ensures a flat rotation curve of the disk due to 
Phaio('') r^^- In this model, the halo dominates L, E and M and the 
disk is a mere tracer of the dynamics. By assuming that the specific 
angular momenta of the halo and of the disk are equal, one can in- 
fer A from an observation of the rotational velocity of the disk, if in 
addition the validity of a Tully-Fisher relation between the disk's 
rotational velocity and mass is assumed, and if a constant factor 
relates the disk mass to the halo mass. Both papers report that the 
measured distribution p(A)dA agrees well with log-normal distribu- 
tion suggested by numerical data, recovering values of (A) ^ 0.05. 
iBerta et al] ( 1200 8t) in addition investigate the influence of spin onto 
the star formation and find an anticorrelation between spin and stel- 
lar mass, as well as the fact that galaxies with recent star formation 
tend to have higher spins. They explain these findings in terms of a 
thinner disk of systems with high values of A, which might be less 
efficient in forming stars, leaving a larger reservoir of gas today. 



5 ANGULAR MOMENTUM CORRELATION 

This section treats the correlation in angular momentum between 
neighbouring haloes: The formalism of the previous section is ex- 
tended for computation of the 2-point variance of the angular mo- 
mentum field and hence the correlation function (Sect. l5.lT l. An im- 
portant simplification, in which the correlations between the inertia 
tensors is disentangled from those of the tidal shear is discussed, 
(Sect. | 5.2|l , follow ed by outlining a very practical model due to 
iLee & PenI ^Ood ) which links the angular momentum directly to 
the tidal shear (Sect. 15. 3t by using a parameter to allow for skew- 
ing between the shear and inertia eigensystems. 



5.1 Ansatz for the angular momentum correlations 

Generalisation of the relations in Sect. 14.11 to include a second 
peak results in the correlation function (L„(x)L„'(x')> of the an- 
gular momenta, with the Gaussian probability density p(w)dw = 
p{v, v')dvdv': 

{L„(x)LAx')) = (66) 
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r duC(w) r &v'C{v') L„(v)LAv')p(v,v'). 

In general, the t hresholds v, v' imposed on the peaks are equal (see 
iLee et aljlioo^ . for a discussion of the additional complication of 
unequal peak heights). As derived in Sect. l3.2l the angular momen- 
tum L„ depends on the product the inertia tensor 7^^^ and the tidal 



shear To- 



(67) 



if the acquisition of angular momentum of a protogalactic object is 
described in the Zel'dovich approximation. Then, the correlation of 
the angular momentum components becomes: 



y{Xpy{x}Xp,y(x'}). 



(68) 



In the next step I replace the Id variance of the components L„ in 
the correlation function by the 3d variance of the full vector L by 
taking the trace of eqn. l |68l ), 

(69) 



Cdr) = liiL„{x)LAx')) = (L„{x)L,{x')), 



which has the advantage of being a coordinate-frame independent 
quantity and allows the usage of the relation e„/3y6, 



Spy'Sp'y for reducing the product of the two e„^y-symbols: 
{L„{x)L„{x')) = a'b\ \(,Xpy{x')Xpy{x)) - <X^/x)X,/,(a:')>] , (70) 

where the order of the indices in the last term is interchanged. In 
matrix notation, the correlation function Cl(?') reads: 

Ci(r) = tr<L(x)L'(x')> = «''i);tr [<X(x')X(x)) - <X(x)X'(x'))] , (71) 

which is non- vanishing for general asymmetric matrices X^y due to 
the matrix transposition in the last term. In the correlation function 
of the angular momenta Cifj) = tr{L' {x)L{x' )) (c.f. eqn. |7U. the 
dependence of L on the commutator X" translates into the asym- 
metric quadratic form (X(jc)X(jc') - X(x)X'(jt:')) with the matrix 
transpose in the second term carrying the signal: A common eigen- 
system of the inertia and shear tensors, being both symmetric ma- 
trices, would have the consequence that X would be a symmetric 
matrix, X = X^ = X', and the correlation function CLir) would 
vanish. 

In contrast, the signal is maximised, if the shear and inertia 
eigensystems are unaligned, i.e. if X is purely antisymmetric, X = 
X" and X^ = 0. In that case X' = (X")' = -X" = -X and the angular 
momentum auto-correlation function 0™'^{r) is truly quadratic: 

C^^''(r) = Id^bl tr [<X(a:')X(x))] , (72) 

which illustrates the effect of partial alignment of the shear and 
inertia eigensystems, reducing the variance compared to the case 
where the shear and inertia eigensystems are maximally mis- 
aligned. Additionally, eqn. J71b yields negative values for antipar- 
allel aligmnen t of the angular momenta of neighbouring haloes (see 
ISchafeJuOOSl . together with a discussion of the symmetry proper- 
ties of eqns.l71landl72t. 



5.2 Decomposition of tidal shear and inertia correlations 

As already pointed out, the correlations in the angular momen- 
tum field result from an interplay of the long-ranged correlations 
in the tidal shear field (which has the same correlation length as the 
density field) and the much shorter range correlation of the inertia 
(which is related to the curvature of the density field, i.e. its sec- 
ond derivative). This means that at large separation distances the 



random processes determining the inertias of the two haloes are 
independent and only take account of the local shape of the den- 
sity field, en coded in the tidal shear. Motivated by this observation, 
ICatelan &~Theuns ( 1996a) propose the decomposition 

(Kxmx) Kx'mx')) = p(\mx)) p{\mx')) {^{x) t(x')), (73) 

with two tensor-valued conditional probability densities 
P{\\'S{x)) = (I ^(jt:)) which give the distribution of obtaining 
a halo with a certain inertia given a particular shape of the density 
field, or equivalently, tidal shear field. This approach is applicable 
on separations larger than the corelation length of the inertia field 
and greatly reduces the complexity of the evaluation of the 4-point 

object <I(x)T(a:) \{x'y¥(x')). 

Treating fo-^ as a continuous field, ICatelan & Porcianil j200lT) 
compute the correlation function ("I'j^yl'o-'y') of the tidal shear field, 
which has to reflect the symmetry relations of the tidal shear tensor 
(by carrying out a decomposition into scalar, vectorial and tensorial 
components), as well as to distinguish between longitudinal and 
transverse correlations with respect to the separation vector. 



5.3 Parameterised angular momentum correlations 

iLee & PenI ( l2000h and lCrittenden et alj ( 1200 ih proposed to link the 
angular momentum directly to the tidal shear and to use an effective 
parameterisation of the misalignment between inertia and shear: 



l+a 



(74) 



where T is the unit-normalised dVap^ap =1) traceless (trT = 0) 
tidal shear tensor, which can easily be derived using 



trT 



•Jap, 



(75) 



and rescaling 4* = 'P/ \^\. The derivation of eqn. MA\ assumes no 
correlation between the inertia and shear tensors, as explained in 
the previous section, uses the decomposition eqn. l l73t applicable 
at large separations and assumes statistical isotropy of the inertia 
tensor, i.e. an evaluation of the inertia correlations in analogy to 
eqn. l [48] l. 



(Ipa-Ip'a->) = (Sf 



+ Spp'Saa' + Sb^'Si 



'Pa-' Op'cr 



+ 4>- 



(76) 



The parameter a, which measures the degree of misalignment be- 
tween the shear and inertia eigensystems assumes values between 
a = corresponding to the case of perfect correlation between in- 
ertia and shear, and the maximum value of a = 3/5, if inertia and 
shear are mutually uncorrelated and if there are no nonlinear con- 
tributions to the angular momentum. If a = 0, one recovers the 
isotropy relation (L„L„/) = S^a'/^, and quite generally, a stronger 
correlation between the shear and inertia makes the angular mo- 
mentum more random. In nu merical simuation s, a was measured 
to assume a value of a = 0.24 jLee & P en '2000"). 

A subtle point is th e parameter isation a in eqn. l l74t : In the 
derivation of the formula, lLee&Penl l l2OO00 assumed uncorrelated 
shear and inertia tensors, as stated above, and proceed by calculat- 
ing the variance of the angular momentum for a given shear, for 
which they find a term quadratic in the tidal shear as a correction 
to the variance resulting from random orientation of the angular 
momenta, because haloes tend to align their angular momentum 
with the intermediate eigenaxis of the shear tensor They argue that 
this dependence on shear should vary between and 1, correspond- 
ing to the cases of random orientation of angular momenta, and of 
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angular momenta fixed by the shear eigenframe, respectively. By 
setting up a Gaussian random process and by calculating the vari- 
ance of the angular momentum direction for a given tidal shear in 
the limit of weak alignment, they recover eqn. J74b . which can be 
scaled to the variance of the full angular momentum by the expec- 
tatio n value (L^ ) /3. 

iLee & Peril j200lh continue by defining their angular momen- 
tum correlation functions r]{r) and rjiir) as the 2-point variance of 
the angular momentum directions L = L/L, 



Ur) = (\L{x)L{x + r)\) - - 



and 



r]{r) = {L{x)L(x + r)), 



1 fl- e^ir) 



(77) 



(78) 



with the filtered correlation function ^R(r) of the density field. 



J 



P{k)W\kR)Ukr). 



(79) 



The proportionality to the square ^\{r) is a direct consequence of 
the parameterisation eqn. l l74t with the squared tidal shear, and the 
term 1/3 makes the correlation function 772(f) vanish in the case of 
no alignment. The length scale R is related to the mass scale M by 
the usual relation M = 4nQ.„,p^y,tR^ /3, applicable for a top-hat filter 
kernel. 

Here, it should be emphasised that eqn. i74l provides a cor- 
relation rather between the angular momentum axes than between 
the angular momenta themselves as it is unable to provide the full 
directional information on the angular momentum vector L and 
can not distinguish between parallel and antiparallel aligned an- 
gular momentum vectors, although it is invariant under a global 
parity transformation. Furthermore, the parameters a and (L~) have 
been measured from collapsed objects in numerical simulations, for 
which the tidal shear mechanism in its simplest form is not appli- 
cable as there are contributions from nonlinear dynamics which ex- 
ceed the description using the ZePdovich approximation. Further- 
more, it is not straightforward in a numerical simulation to iden- 
tify and to measure the initial conditions from which an object has 
formed. Additionally, eqn. i74l links the correlation properties of 
the continuous tidal shear field to those of the discrete angular mo- 
mentum field, without the biasing of angular momentum. Never- 
theless, the formula has enormous merits in describing correlations 
in the angular momentum direction, and is successfully used in in- 
trinsic shape alignment studies, as the galaxy ellipticity does not 
depend on the magnitude of L and is invariant under the transfor- 
mation L — > -L. 



6 NUMERICAL SIMULATIONS 

Angular momentum build-up by tidal torquing and the result- 
ing angular momentum correlations have been extensively tested 
in numerical «-body simulations of cosmic structure formation 
(Sect. 16. U , as well as alignments of angular momenta with am- 
bient cosmological structures such as filaments, sheets and voids 
(Sect. |6.2t . In addition, there are studies that focus on compet- 
ing mechanisms for angular momentum build-up such as accretion 
(Sect.l63t. 



6.1 Numerical tests of tidal torque theory 

One of the first numerical studies of angular momentum build- 
up of objects in t he lar ge-scale structure was carried out by 
Efstathiou & JonesI j 19791) . who ran a simulation of 10' particles, 
finding a few tens of objects. They were able to reproduce the L oc 
M^'^-relationship predicted by linear tidal torquing, and analysed 
the distributions of the spin parameter A, reco vering the low values 
for A p redicted analytically. Later papers, e.g. lBamer & Efstathio^ 
improve conclusions tremendously due to the increase in 
number of particles, realistic initial conditions and algorithmic ad- 
vancement. Their results include an accurate determination of the 
L-M-relation, where they find best fit values close to the fiducial 
L oc M^'^-scaling, as well as the time evolution of the angular mo- 
mentum acquisition, with its truncation at turn-around. The distri- 
butions of A were found to be insensitive on the particular cosmol- 
ogy and on the slope of the matter power spectrum P(k) oc k"' . 
Additionally, they investigated the alignment of the spin axis with 
the ellipticity axes of the dark matter halo and found the spin to 
point preferentally perpendicular relative to the major axis of the el- 
lipsoid, before they discuss the relation between A and mo rphologi- 
cal typ e, emphasising the importance of baryonic phvsics. IPrimackl 



J— ^' - - - 

12003) provides a summary on the relations between halo angular 

momentum and disk formation, emphasising open questions on the 

matter. 

ISugermanetalJ ( l2OO0l) carry out an extensive comparison be- 
tween the analytical predictions from tidal torquing with numerical 
n-body simulations, confirming the basic tidal torquing mechanism. 
They find a linear growth with cosmic time until turn-around, and 
continues quasi-linearly until shell crossing, when the Zel'dovich- 
approximation ceases to be applicable. Linear tidal torquing was 
found to overestimate the true turn around by a factor of three 
(indicating the importance of tidal forces on object formation by 
collapse) as well as yielding too large values for the total angular 
momentum. In measuring the mass and binding energy of the dark 
matter haloes, they confirm the anti-correlation of the spin param- 
eter A with peak height v and the angular momentum scaling with 
mass, L oc M''\ 

Details of tidal torquing as a mechanism of angular momen- 
tum build-up are scr utinis ed in the seri e s of pa pers written by 
IPorciani et all ( l2002al) and IPorciani et al.l j2002bl) . In the first pa- 
per, they compare predictions from tidal torquing to the spin am- 
plitudes of haloes forming in a cosmological ;7-body simulation, 
assuming that torquing is efficient up to 0.6 halo's turn around time, 
and find that the order of magnitude of the angular momentum can 
be well predicted down to redshifts below z - 3- Nonlinear influ- 
ences cause a significant deviation of the angular momentum direc- 
tion from the direction predicted by tidal torquing, amounting to a 
mean error of about 50°, which weakens the spin correlations, and 
casts doubts on the usefulness of tidal torquing for pre dicting ellip- 
ticity alignments. In particular. jPorciani et al. ( l2002al) confirm that 
the mechanism of tidal torquing in the linear regime is a very good 
model, and that the truely nonlinear effects at late times cause de- 
viations: Not only is the rotation axis affected, but the amplitude of 
the angular momentum grows by more than a factor of two between 
redshifts of z = 3 and z = in addition to the prediction of linear 
tidal torquing. Further results include the angular momentum time 
evolution in units of the turn-around time to, which levels out from 
a linear increase at =^ 0.6tQ, in accordance with the tidal torquing 
model, as well as the distribution of ellipticity e and prolateness 
p, which agrees well with the expected distribution for a Gaussian 
random field. 
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Figure 4. Alignments of the tidal torque and inertia eigentsystems, quanti- 
fied by the distributions of the cosines of the angles between the principal 
axes of the tidal shear and in ertia tensors. T h e figure was kindly provided by 
C. Porciani and appears in \Porciani et al\ k002ti) . Reproduced with per- 
mission from Blackwell Publishing. 



Two popular ways of defining spin correlation functions of the 
angular momentum directions L = L/L are 

ri(r) = {Ux)Ux + r)), (80) 

ri2(r) = {\L(x)L{x + r)\) - ^, (81) 

which have been m easured from the n-body simulation by 
iPorciani et al. who find a positive spin correlation func- 

tion at early times out to separations of a few comoving Mpc, but 
this correlation decreases towards z = 0, due to nonlinear effects, as 
shown in Fig. |3] In their second paper on the topic, Porciani et al., 
( l2002bl) focus on misalignments between the tidal shear and inertia 
tensors, and investigate the extend to which the common assump- 
tion of uncorrelatedness is fulfilled. To this end, they determine the 
eigenvectors and /„ of the two tensors and measure the distribu- 
tion of the cosines between two pairs of vectors. They find strong 
alignment between the minor axes of the two tensors, as illustrated 
in Fig. |4] which gives the distributions of the angles between ^„ 
and /g. The correlation parameter for the directional cosines yu is 
found to be ^ 0.61, again indicating strong correlation between 
shear and inertia. This correlation would affect the decomposition 
eqn. l [73] l. 

A test of the models for nonlinear torquing and the result- 
ing spin correlation functions from numerical data is provided 
by [L ee & Pen (2007b). Their starting point is the observation by 
iHui & Zhana (2002), that the correlation r]2{r) oc ^-(r) is only valid 
for Gaussian fields, and introduce a correction oc f(r) to the cor- 
relation function, which implies that the spin correlation is more 
long-ranged if nonlinear modifications to the dynamics are present. 
In fact, their data shows weak correlations out to 20 Mpc//i, com- 
pared to 10 Mpc//j for linear tidal torquing, with nonlinear effects 
increasing with decreasing redshift, and is well fitted by the model 
with the additional correction term. 



6.2 Angular momentum alignments with the cosmic web 

Apart from the angular momentum correlation function, which is 
due to correlations in the initial conditions for angular momentum 
acquisition, there are alignments between angular momenta and the 
ambient large-scale structure caused by the large-scale tidal gravi- 
tational fields o f nonlinear structure s like sheets and filaments (for 
an overview, see lColberg et al.l2005h . which has attracted consider- 
able interest. Due to domination of nonlinear physics, these align- 
ments can only be investigated in n-body simulations. Primary top- 
ics of the analyses are, apart from angular momentum alignments, 
the alignments of the halo axes with the large-scale structure, and 
the orientation of galactic disks. 

[Nav arro et al .l |2004|) focus on the orientation of the angular 
momentum axis of disk galaxies (i.e. the angular momentum of 
the dark matter haloes of disk galaxies) and confirm the expecta- 
tion that the haloe revolves around the intermediate axis. Extend- 
ing their analysis they investigate the orientations of the galactic 
disk and of the halo relative to the surrounding large-scale struc- 
ture, they find that the orientation of rotation axis lies in a plane 
traced by the protogalactic material, which provides an explanation 
of the high inclination of Milky Way relative to the supergalactic 
plane. 

The paper bv lSailin & Steinmetj ( |20()^ addresses the align- 
ment of shapes and anngular momenta inside a halo as well as rel- 
ative to the ambient density field. Their results on the first point 
include the an alignment of the angular momentum vector with 
the intermediate axis of the halo ellipsoid, which shows a slight 
dependence on radius. The external alignment of the halo's axes 
and of the angular momentum was quantified by correlations be- 
tween the principal axis system and the separation vector between 
neighbouring haloes, for which evidence was found, correlations 
between the principal axis systems of two haloes, which suggests 
correlations between the major axes only, correlations between the 
angular momenta and the separation vector between two haloes, 
with no conclusive evidence and finally the angular momentum 
correl ation funct i on, wh ere no deviations from isotropy were mea- 
sured. iLee et al] l l2007h choose to quantify shape correlations by 
ellipticity correlation function and ellipticity-direction correlations, 
and carry out measurements of those correlation functions, finding 
strong support on tidal torquing, and providing polynomials in the 
density correlation function as fitting formulae. 

IPatiri et alj ( l2006h and iBrunino et al] ilOOlh investigate the 
alignment of dark matter haloes with voids in the large-scale struc- 
ture and provide evidence that the minor axis of dark matter haloes 
points preferentially to the centres of voids, while the angular mo- 
mentum was not found to have a particular orientation. Conse- 
quently, the intermediate and major axes were found to point pref- 
erentially perpendicularly to the void direction. These results ap- 
ply to voids larger than 10 Mpc /A for haloes in a shell of thick- 
ness 4 Mpc/h around the void. The directional cosines yu be- 
tween the void direction vector and the principal axes of the el- 
lipsoid can be described by a probability distribution of the shape 
P(ji) oc p/(l -I- (p^ - l)yu)^'-, where p = 1 corresponds to a case of 
no alignment. 

The topic of the analysis bv lAragon-Calvo et al.l l l2007h is the 
spin alignment with nonlinear features in the large-scale structure. 
They use a filtering scheme for segmenting the large-scale structure 
into sheets and filaments, using the curvature tensor djdjd of the 
density field, which they reconstruct from the particle distribution 
with a Delauney tesselation technique. The ratio of eigenvalues of 
djdjS provide measures of local spherical symmetry, filamentarity 
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Figure 3. Correlation function and Tj2{r) of the angular momentum directions of low high mass (left) and low mass (right) haloes, comparing the predictions 
from tidal torquing (open squares) with the results from ra-body simulations (filled circles). The error bars reflect sampling errors only and ignore the error 
associated wit h the spin measurement. The dashed line indicates the theoretical expectation. The figure was kindly provided by C. Porciani and appears in 
\Porciani et al\ ^2002bJ. Reproduced with permission from Blackwell Publishing. 



or planarity: Haloes in walls have spin vectors lying in the plane of 
the sheet they are situated in, the orientation of spins in filaments is 
weaker, but acquires a mass dependence. 

'Hah n et alj ( l2007h investigate ellipticity and spin-parameter 
distributions, as well as the spin-spin and spin-orbit correlation 
functions, studying the enviroimient and redshift dependence of 
these observables. As environmental classification they use the 
number of positive eignvalues of the (smoothed) tidal shear ten- 
sor, which is motivated by the fact that compression of the cosmic 
density field takes place along directions where the tidal shear has 
a positive eigenvalue. For the eigenvalues /I,, ; = 1, 2, 3 of the tidal 
shear one can distinguish between clusters (/ii,2,3 > 0), fil- 

aments (/i2,3 > Q, A[ < 0), sheets {X^ > 0, A12 < 0) and voids 
Wi,2,3 < 0). Their results comprise trends of halo sphericity and 
triaxiality with environmental classification, as well as enviroment- 
dependencies of the spin parameter A, which on average assumes 
higher values in filaments for high mass haloes, and significantly 
lower values for low mass haloes in voids. 

Spins of low-mass haloes tend to be aligned with the filament 
direction, haloes similar to the nonlinear mass scale are randomly 
aligned whereas high-mass haloes have a weak tendency to orient 
their spin perpendicularly to the filament. ^ Hahn et aL 1/2007) did not 
find a spin-spin correlation between neighbouring haloes, but re- 
port a substantial spin-orbit correlation, independent of mass-scale 
or environment: Pairs of haloes tend to align their spins parallely 
to the orbital angular momentum. Fig. [5] shows the tidal fields and 
their alignment relative to the large-scale de nsity field in an n-body 
simulation carried out bv lHahn et al.l JioOTl) . indicating the impor- 
tant influence of nonlinear dynamics on the tidal shear. 

ICuesta et al.l l l2008l) revisit the spin alignment with voids. 
They confirm the v alidity of the fitting formula proposed by 
iBrunino et al.l ( |2007|) for the angle distribution between the ellip- 
soid's principal axes and the void direction ruling out random align- 
ment with high significance, and in addition, find a preferential 
alignment between the angular momentum and the void direction. 




h 'Mpc 



Figure 5. Slice through a simulation of the large scale density field: Grey 
dots indicate haloes, and the arrows denote the direction of the eigenvec- 
tor corresponding to the largest eigenvalue of the tidal shear tensor inside 
filamentary regions, and is used for measuring shape or angular momen- 
tum alignment. The figure was kindly provided by O. Hahn and appears in 
[Hahn et al. QOOTJ. Reproduced with permission from Blackwell Publishing. 



which was not present in earlier studies. They attribute this result to 
the criteria by which particles were asssigned to haloes, which im- 
pacts on the determination of the inertia tensor and of the angular 
momentum. 
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6.3 Spin-up from merging and satellite accretion 

Apart from tidal torquing a galaxy may acquire angular momen- 
tum by accretion of satellites. These satellites may bring in their 
own angular momentum or enter the merging system under a cer- 
tain impact parameter, such that in the process of accretion orbital 
angular momentum of the satellite gets transformed into spin an- 
gular momentum of the merged system. The spin-up in s atellit e ac- 
cretion was considered among others bv lVitvitskaetalJj2002h and 
iD'Onghia & Navarrol ( l2007h who find in n-body simulations that 
the mechanism of sa tellite accretion is sufficient to provide enough 
angular momentum. iD'Onghia & Navarr design a model 

for angular momentum build-up my merging, which is based on 
stoachastic increase in angular momentum in the sequence of merg- 
ing processes of a dark matter halo, and which is able to reproduce 
the log-normal distribut ion of the spin parameter A . Extending on 
the preceeding work bv lD'Onghia & Burkerj ( l2004h who find that 
tidal torquing does not provide enough angular momentum for ro- 
tationally supported objects to f orm and conclude that accretion as 
a further mechanism is required. lD'Onghia & Navarro ( 2007) show 
that the distribution of A depends on the dynamical state of the halo 
and shows difl^erences between relaxed haloes and haloes which 
underwent recent merging. The relaxation of a halo is quantified by 
the distance i- between the most bound particle and the centre of 
mass, which should assume small values for a relaxed, spherically 
symmetric object. They quantify the time evolution of s, A and halo 
mass^ 

iPeirani et al. Il l200i) take these investigations further and clas- 
sify the systems according to their merging history into two cate- 
gories: continuous accretion and sigificant merging. Their results 
suggests that the angular momentum growth exhibits a transition 
time at z 1.5, where the fast initial growth levels off. In addition, 
they find a correlation between the mass and angular momentum 
evolution, suggesting that merging has a significant impact on the 
angular momentum build-up history. The spin parameter A is de- 
creasing with time for haloes undergoing accretion, which is com- 
pensated by a small increase in merging haloes. 

The numerical studi es are supplemented b y an analytical in- 
vestigation carried out bv lOuinn & BinnevI ( IT992) . who considered 
angular momentum re-orientation in an accretion model. They con- 
struct spherical shells of matter, average the torque of each of those 
shells as a function of the density and potential, and express those 
two fields in terms of multipole moments, finding that the multi- 
pole series converges quickly such that only the di- and quadrupole 
moments contribute. They conclude that especially the dipole term, 
which gives rise to anticorrelated torques in neighbouring shells, 
and gives rise to warps as the galaxy can not maintain a common 
rotation axis throughout the dark matter halo. 



7 OBSERVATIONS 

Intrinsic ellipticity correlations were measured using various statis- 
tics in a number of cosmological surveys. For instance. Pen et al. 
( l2000h measured the spin-direction correlation (|L(xi) • r(x2)| ), 
r = X2 - Xi, in the Tully galaxy catalogue, an all-sky survey com- 
prising 3.5 X 10'' galaxies out to redshifts of z 0.02. Although 
the signal is weak, the measured correlation strength is consis- 
tent with theoretical estimates. The spin-spin correlation function 
<|L(xi) ■ Lfe)^ > is measured bv lLee & PenlfeOOll) on a subsample 
of the same data set, yielding a total significance of 2Aa. They 



confirm that the simplified angular momentum coupling model 
eqn. J74[> with the param eter a ^ 0.24 fitting the data nicely. 

Brow n et"al] ^200 2^ set out to measure the alignment of ellip- 
tical galaxies in the SuperCOSMOS data set, which reaches out 
to a median redshift Znied = 0.1 and comprises 2 x 10* galax- 
ies over a quarter of the sky. As a statistical quantifyer, they use 
the variance of the obseved ellipticity, and show that it is larger 
than the expected lensing signal, and well in accordance with the 
predicition from quad ratic alignment models, in particular with 
ICrittenden eTal] ( 1200 ih . 

Another high-sig nificance measure ment of the spin-shear cor- 
relation is reported bv lLee&Pej ( l2002h . who constructed both the 
spin field as well as the tidal shear field (via the density field) from 
the PSCz-catalogue, a complete redshift catalogue from the IRAS- 
survey covering 84% of the sky, and directly quantified the correla- 
tion between the two fields. The null-hypothesis of random align- 
ments is rejected at 99.98% confidence, and the alignment parame- 
ter was found to be a 0.17 ± 0.04, slightly smaller than the value 
suggested by numerical simulations, although the authors emphasis 
that their data treatment i s likely to systematically underestimate a. 

iTruiillo et al.l j2006h perform a measurement of the orienta- 
tion of the rotation axes of galaxies situated on the shells of large 
voids, using the 2dFGRS and SDSS data sets. They measure the 
distribution of the angle between spin axis and the density gradient 
(i.e. a related quantity to the void centre direction) and can reject 
the null hypothesis of randomly aligned galaxies at 99.7% confi- 
dence. Comparing their distribution to the prediction derived from 
eqn. MA\ , the data requires much larger values for a, a ^ 0.7, but 
prov ides otherwise an excel lent fit. 

iLee & Erdogdul j2007h combine a reconstruction of the den- 
sity respective tidal shear field from the 2MASS reshift survey 
with the spin field derived from the Tully galaxy catalogue. The 
measured correlation strength exceeds 6cr significance, and allow 
the investigation of the morphology and environment dependence. 
iLee & Erdogdul ( 12007 ) find a weak increasing trend of the corre- 
lation with increasing morphological type and measure larger cor- 
relations in high-density environments. Additionally, they confirm 
the predicitions of tidal torquing to be well describing the data, in 
particular the probability desity distribution of the cosine between 
the spin and the intermediate axis of the tidal shear tensor. Their 
mean correlation parameter a =^ 0.084 ± 0.014 is much smaller than 
the value obtained in numerical simulations. 

Aiming at distinguishi ng different alignm ent scenarios for el- 
liptical and spiral galaxies. ILee & PenI ( l200 7al carry out an inves- 
tigation of alignments of red and blue galaxies in SDSS, selecting 
galaxies out to redshifts of z = 0.4. For the description of intrinsic 
shape correlation the functional form ri(r) oc a^^(r) + €^{r) is used, 
where the case a <c 6 would correspond to the usual quadratic 
alignment model for spiral galaxies with a small correction due 
to non-Gaussianity, whereas a » e would be applied to ellipti- 
cal galaxies, where a linear alignment model applies. The measure- 
ment in fact supports this ideas as it yields a = 0.20 ± 0.04 and 
e ~ for blue galaxies and a - and e = (2.5 ± 0.5) x 10"' for red 
galaxies. 



8 IMPLICATIONS FOR GRAVITATIONAL LENSING 

An important application of angular momentum correlation models 
is the problem of intrinsic alignments in weak gravitational lens- 
ing. After a short recapitulation of gravitational lensing (Sect. lSTt 
and especially weak cosmic shear (Sect. 18. 2t . I explain the con- 
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tamination of lensing data by intrinsic aligments in tiie measure- 
ment of weak sfiear spectra (Sect. [83] ) and review common models 
linking ellipticity correlations to angular momentum correlations 
(Sect. l8.4] l. Then, I discuss measurements of the intrinsic alignment 
terms, both on simulations as well as on actual data (Sect. [83] l and 
give an overview over methods proposed to suppressing intrinsic 
alignments (Sect. lS^ . 

8.1 Gravitational lensing 

Gravi tational lensing (for a review, see 'Bartelman n & Schneidej 
l200lh refers to the deflection of light reaching the observer from 
a distanct galaxy by intervening gravitational potentials. The posi- 
tion /3 at which a galaxy is observed is related to the angle under 
which the galaxy would be observed without gravitational lensing 
by the lens equation. 



(82) 



with the deflection field a(0), which is given as the angular gradi- 
ent a = Vif/ of the deflection potential if/, if/ follows from the surface 
mass density 1.(9) of the lens as a solution to the 2d Poisson equa- 
tion. 



Atf/(e) = 2k(9) 11/(9) ■■ 



J d-0' Ki 



(6»')ln|6»-6»'|, 



(83) 



after normalisation with the surface critical density Scrit = 
l(AnG)Dsl(DiiDijs) such that the convergence k is related to the 
surface mass density of the lens via k(9) = E(0)/E(.,it. The angular 
diameter distance to the lensed galaxy is denoted by and like- 
wise Dd is the angular diameter distance to the lens and D^, the 
distance between lens and background galaxy. Since 9 can princi- 
pally never be known, one rather observes linear variations of the 
deflection field a across the image of the galaxy leading to shape 
changes in the galaxy image. This information is contained in the 
Jacobian 



89' 



deidOj 



(84) 



which describes the shape distorsions of an image to first order. The 
difference A,?l between the Jacobian and the unit-mapping can be 
decomposed with the Pauli-matrices erg, 



11=0 



(85) 



because the cr^ constitut e a complete basis of th e vector space of 
the real 2 x 2-matrices jArfken & Webeill2005h . The coefficients 
are refered to as the convergence k, which describes the isotropic 
image size change, the two components of shear y+ and yx, and 
the image rotation p, which can only be excited by multiple lensing 
events along the line of sight or if corrections due to the perturbed 
photon geodesic become appreciable. Due to the property cr„(rp = 
id(2)dafi + ie„^yO"-y of the Pauli-matrices, the coefficients can be re- 
covered using cip = tr(A^cr^)/2. Due to their common transforma- 
tion under coordinate system rotations, the two shear components 
can be combined to the complex ellipticity 7 = 7++ i7x, which 
forms a spin-2 field and transform according to 7 — > 7 exp(2i(fi) in 
rotations about an angle ip. 

The shape of a galaxy is measured from the second moments 
of the brightness distribution Qjj, which in turn defines the complex 
ellipticity e, 

211-622 , 22,2 



with 



Q 



and 



Q 



(86) 



where I abbreviated 2 = 2ii + 222 + 2(2ii 222 - 2^2) The com- 
plex ellipticity forms a spin-2 tensor field due to the transformation 
property 6 — > eexp(2i(^) under the rotation of the coordinate sys- 
tem by an angle (p and maps onto itself if (f = n. Under the action 
of gravitational lensing, the intrinsic ellipticit y 6j of a galaxy is 
distor ted to the observed ellipticity 6 given by jSeitz & Schneiden 
[1993) 



l+8'es 



(87) 



with the complex reduced shear g = -//(l - k). In weak lensing 
applications the relation simplifies to g ^ 7 because \k\ <g: 1. 



8.2 Weak cosmic shear 

Weak gravitational lensing directly probes the fluctuations of the 
cosmic density field by imprinting correlated changes to the shapes 
of neighbouring galaxies, because the light rays reaching us from 
pairs of galaxies have to transverse the same part of the large-scale 
structure and experience correlated tidal fields. In this context, the 
term weak refers to the fact that the distorsions imprinted onto the 
galaxy shapes are small, k <g: 1, I7I <c 1, and the signal is obtained 
as a statistical average over many galaxies. Although the primary 
observable is the shearing of the galaxy image, one often works 
in terms of the lensing convergence k as it has the same statistical 
properties as the shear and is easier to work with. 

The convergence k is given by a line of sight integration: 



AxG(x)x—5, 
Jo a 



(88) 



replacing the actual photon geodesic by a straight line (Bom- 
approximation), with the lensing-efficiency weighted distribution 
of galaxies G(x), 



G(x) 



r 



dy p(x') 



x' 



(89) 



The redshift distribution of the lensing galaxies is denoted by 
p(x)dx, after changing the distance measure from redshift to co- 
moving distance. The angular power spectrum Q(^) of the lensing 
observable is related to the power spe ctrum P(k) of the density field 
via Limber-projection ( lLimbeill954) . 
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where the integration range formally extends to the comoving hori- 
zon distance j^'/^. Technically, one can tremendously improve the 
sensitivity of weak lensing by dividing the galaxy sample into red- 
shift bins and carry out weak lensing tomography, where the lens- 
ing spectrum is measured separately from each redshift bin con- 
taining the respective information on structure growth and the geo- 
metrical factors, and by combining those measurements. The weak 
lensing convergence /c, in a tomographic measurement reads 



D+ 

Ki = dxGi(x)x — 
Jo a 

with the altered lensing efficiency function G,(^) 

G,{x)= r dx' Pi(x')^^-^, 

Jx X' 



(91) 



(92) 



constructed for the galaxy distance distribution Pi(x)dx of the red- 
shift bin i. In this way, on can access cosmological information 
that would be otherwise diluted by the line of sight integration 
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eqn. ( I88t with a slowly varying w eighting function. A f urther de- 
velopment of this idea proposed bv lHeavens et aL H2006h is refered 
to as 3d-lensing, where the density field and its time evolution is 
reconstructed from the weak lensing shear, which is effectively a 
measurement of the local tidal field. 



8.3 Intrinsic alignments 

A common assumption in weak lensing are uncorrelated intrinsic 
shapes, which is challenged by the fact that there are ellipticity 
correlations as a consequence of angular momentum correlations. 
In the case of uncorrelated intrinsic ellipticities, the observed an- 
gular convergence power spectrum differs from the theoreti- 
cally expected spectrum C,(f) by a Poissian noise term cr^/n, which 
eventually dominates on small angular scales, 

2 

CM) = €,{() + —■ (93) 
n 

The dispersion of the Gaussian ellipticity distribution is usually 
quoted to be cr^ = 0.4, and n denotes the number of galaxies per 
steradian. Typical number densities of galaxies range between 20 
and 100 objects per squared arcminute, depending on the survey 
setup. If the galaxy shapes are intrinsically correlated, the above 
relation becomes 



CM = CM + CM- 



(94) 



with the intrinsic ellipticity correlation Ci(£) if the intrinsic shapes 
are uncorrelated with the lensing shear. This is actually not the case, 
because the same tidal shear field gives rise to both the intrinsic 
ellipticity alignment as well as the lensing signal, i.e. there can b e 
an additional cross correlation term Qe(^) ( iHirata & Seliakll2004l) . 



CM = CM + CM) + C,M, 



(95) 



especially in linear alignment models. To be exact, Cji^l) describes 
the correlation between the lensing induced apparent shape change 
of a background galaxy and the physical alignment of a foreground 
galaxy in the tidal field of the gravitational lens. The physical shape 
of a background galaxy can safely be assumed to be uncorrelated 
with the tidal field of the lens as they are typically separated my 
many correlation lengths. 

The quantification of intrinsic ellipticity correlations in SDSS 
an d the implications for wea k lensing is the topic of the paper 
by iMandelbaum et all | |2006|) : Using a spectroscopic sample of 
2.6 X 10' galaxies at low redshift (z i 0.12), they carry out mea- 
suremets of the intrinsic ellipticity correlation CM) as well as of 
the cross correlation C^M) between ellipticity and density. While 
the first correlation was too week to be detected, the latter correla- 
tion, present in linear alignment models, yielded a clear detection 
in bright galaxies out to scales of 60 Mpc//i, with the correct sign. 
Inferring consequences for weak lensing from their results suggests 
that the contamination of deep lensing surveys with intrinsic align- 
ments could be as large as 10%, and that values of erg derived from 
weak lensing might be overestimated by as much as 20%, if un- 



corrected. Addition ally, they provide fits to the alignment model 

lOOQ)forcar i 

Extending this analysis, Hirata et al] j2007l) investigate corre- 



, they pri 

proposed by Heavens et al.l j200Q) for carrying o ut that correction. 



lations between ellipticity and the density field, using 3.6 x 10' 
LRGs from SDSS and 8000 LRGs from 2SLAQ, resulting in a > 
3cr detection of the alignment effect on large scales up to 60 Mpc//;. 
The high detection significance allows the determination of fitting 
functions for the scaling of alignments in relation to luminosity, 
transverse separation and redshift. Applying their findings to the 



estimation of cosmological parameters from weak lensing surveys 
they conclude that the contamination can decrease the lensing spec- 
trum by typically 6.5%, which induces as bias in the estimation of 

erg o f Ao-g = -0.02 

iBridle & Kind ( |2007|) show that the intrinsic alignment con- 
tributions to weak lensing spectra cause severe errors in the esti- 
mation of cosmological parameters, in particular the dark energy 
equation of state parameter w, which can b e biased by as much as 
50% if the alignments remain uncorrected. iBridle & K in'3 ll2007h 
emphasise the importance of accurate photometry in tomographic 
weak shear measurements, with respect to the total number of red- 
shifit bins as well as the photometric redshift accuracy. 



8.4 Ellipticity correlations 

A consequence of angular momentum correlations between neigh- 
bouring spiral galaxies are the correlations between their apparent 
ellipticities, because their galactic disks are viewed under similar 
angles of inclination. In doing that, one assums that the galactic 
disks form perpendicular to the angular momentum axis of the host 
halo. For disks with zero thickness the relation between the scalar 
ellipticity e and the line of sight component of the angular mo- 
mentum is given by 

1-L? 

6= TT- (96) 

\+Ll 

Being a purely geometrical effect, the ellipticity depends only on 
the direction L = L/L and not on the absolute value of the angular 
momentum. The two components of the complex ellipticity e = 
6+ -I- i6x can be derived from the angular momentum direction via 

?2 _ f2 
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and 
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such that e- = 6; + ej. This alignment model is commonly refered 
to as the quadratic alignment model, due to the dependence on Lr. 
Contrarily, the alignment model applying to elliptical galaxies is 
usually parameterised to be linear in L. 

A finite thickness of the galactic disc can be incorporated by a 
constant of proportionality a in the range < a < 1, 



' 1 + L? ■ 



(98) 



which weakens the dependence of e on L-. ICrittenden et aklfeOOlh 
measure a to have a value of a- ^ 0.85 for spirals, and a - 0.5 
for spheroids and ellipticals. For typical mixture between spiral, 
lenticular and elliptical galaxies from the APM-survey, they obtain 
a mean value of a ^ 0.73. 

After li nking the ellipticity corr elations to angular momentum 
correlations. ICrittenden et aU ( 1200 11) continue by using the decom- 
position eqn. i73\ along with multivariate Gaussians for the condi- 
tional probability density P(L\W) and the quadratic model eqn. i74l 
for the relation between angular momentum and shear for deriving 
the 2-point correlation function of the ellipticity field. They con- 
firm the important result that the ellipticity correlation function is 
proportional to the squared correlation function ^(r) of the density 
field. In the next step they use a Limber-projection with a weight- 
ing proportional to the galaxy density including galaxy clustering 
in order to obtain the angular correlation function of the intrinsic 
ellipticity, which then can be compared to the angular correlation 
function of the weak lensing shear. Intrinsic alignments were found 
to be a serious contaminant only for shallow surveys, as illustrated 
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in Fig.|6l A survey with median redshift of 0. 1 would be dominated 
by the intrinsic alignment signal, whereas the intrinsic ellipticity 
alignments are smaller by an order of magnit ude for a deep sur- 
vey with a median survey redshift of unity. Le e & Peril | |200 (J) take 
this result to extremes and show that the density field could be re- 
constructed from the ellipticity field due to the angular momentum 
distribution instead of using the tidal field measured by weak lens- 
ing, by deriving a formul a for the spin-d ensity correlation. 

In a related work, iMackev et ai] ^2002) investigate angu- 
lar momenutum-ind uced ellipticity correlat ions, using the same 
quadratic model as ICrittenden et alj ( 1200 derive the angular 
power spectra and compare the signal amplitude to the lensing 
power spectra. Technically, they treat the tidal field as being con- 
tinuous and the inertia and shear as uncorrelated. Comparing the 
intrinsic alignment to the cosmological lensing signal, they con- 
firm the dominance of the intrinsic alignments over the lensing for 
surveys with low source redshifts of z ^ 0.1, equality of both ef- 
fects for redshifts of z ^ 0.3 and show that a deep survey with z- I 
has a negligible contamination of the lensing signal by more than 
two orders of magnitude. 

In what concerns galaxy-galaxy lensing, where the foreground 
galaxy density is correlated against the background ellipticity field 
at small angular separat ions, the impact of in trinsic alignments is 
investigated in detail bv lHui & Zhand j2002h . Their starting point 
is the acute observation that tidal models, where the alignment 
is proportional to the squared tidal shear, can not cause a non- 
vanishing ellipticity-galaxy density cross-correlation {ye): The cor- 
relator could be expanded into a 3-point function, which would 
vanish, provided that the density field is Gaussian, the bias is unity 
and there would be no source clustering. They continue by quanti- 
fying the magnitude of the correlation (ye) for non-Gaussianities 
which arise in the course of structure formation. While the re- 
sult is very interesting for cosmological weak lensing and for test- 
ing intrinsic alignment models in themselves, they bridge the gap 
to galaxy-galaxy lensing by showing that there are influences by 
intrinsic alignments because the separation into fore- and back- 
ground objects is imperfect. An investigation of the transition be- 
tween galaxy-galaxy len sing and cosmic weak shear is provided by 
iBridle & Abd5i3 bOOTh . 



8.5 Quantification of the intrinsic alignment signal 



One of the first studies for quantifying intrinsic alignment in com- 
parison to the cosmological weak lensing signal was carried out by 
iHeavens et al.l (12000.), who measured the ellipticity spectrum C^it) 
from haloes in an w-body simulation, assuming that galactic disks 
are aligned with the angular m omentum of the host halo. Quite in- 
dependent of cosmology. Heav ens et alj ( |2000|) found that the weak 
shear signal Ci^it) dominates over Q (f) by about an order of mag- 
nitude, for a deep survey, but show that the reverse applies to shal- 
low surveys, for angular momentum aligned spiral galaxies. For the 
case of elliptical galaxies, which are assumed to have the same axis 
system as the host halo, they report the measurement of an ellip- 
ticity correlation function with an amplitude twice as large as for 
spiral galaxies, but show that the lensing signal of a deep survey 

dominates over this intrinsic elliptici ty correlation. 

Using a very different model, ICatelan et al.l ( 1200 ll) provide 
a theoretical discussion of intrinsic alignments in comparison to 
the weak lening signal, by assuming that the galaxies are distorted 
physically by the tidal shear field, which causes the ellipticity to 
be proportional to the tidal shear and not the square. They recover 
the result that deep surveys are almost free from intrinsic alignment 



contamination, as the spectra differ by almost two orders of magni- 
tude, and that intrinsic alignments become appreciable in shallow 
surveys. Comparing their alignment model with alignments due 
to angular momentum correlations they emphasise that angular- 
momentum induced ellipticity alignments, being proportional to 
the square of the tidal shear, are short-ranged compared to align- 
ments directly caused by ellipticity deformations due to the tidal 

shears 

ICroft&MetzledJ2000l) measure the correlation of ellipticities 
of dark matter haloes in an n-body simulation, assuming that the 
halo shapes serve as proxies for the galaxy shapes, which are mea- 
sured from the second moments of the projected mass distribution. 
They report a positive correlation out to distances of 20 Mpc/A, 
and find that the signal is not strongly affected by varying the 
halo-finding technique (in particular the friends-of-friends linking 
length), nor by the numerical resolution. In comparing the intrinsic 
ellipticity correlation to the weak lensing signal, they show that 
the ellipticity correlation is a small contamination for deep sur- 
veys, but mimicks the functional behaviour of the lensing signal. 
Additionally, they confirm that the first reports on measurements 
of cosmic shear actually pick up the cosmological weak lensing 
signal, and that the intrinsic ali gnment contaminat i on amounts to 
rough tly 10% in these surveys i Bacon et al. 20001 : Iwittman et al.l 

20001). 



120001 : iKaiser et alj200(]| : IVan Waerbeke et al 



2002h developed in- 



Again working with n-body data, iJing 
trinsic alignment studies for weak lensing applications: The ob- 
served ellipticity was assumed to be determined by the halo ellip- 
ticity, which allows the direct measurement of the ellipticity corre- 
lation function as well as of the angular ellipticity correlation func- 
tion. Placing the haloes at unit redshift a nd comparin g with the 
lensing convergence correlation function, Ijind ( 12002? ) confirmed 
that intrinsic alignments may be a serious weak lensing contami- 
nant. It should be emphasised, however, that in these early stud- 
ies an alignment model applicable to elliptical galaxies was used, 
which was later shown to generate a significantly higher correlation 
signal compared to a model based on angular momentum align- 
ments, which is a better choice for spiral galaxies, because they 
make up the majority of field galaxies. In particular, an alignment 
model based purely on halo-ellipticity would be already ruled out 
by weak lensing o bservations, as it overp redicts the signal. 

The work bv lHirata & SeliakI j2004h compares the magnitude 
of intrinsic ellipticity correlations and ellipticity-shear cross corre- 
lations in a simplified linear alignment model, and investigate the 
dependence on source redshift. The usage of the linear alignment 
model is motivated by the fact that quadratic models do predict 
the cross-correlation between ellipticity and shear to vanish: In the 
framework of Gaussian fluctuations, one would need to compute 
the expectation value of an odd moment of the density field, which 
vanishes exactly. In an analytical calculation they derive that in the 
linear alignment model the primary contaminant of the weak lens- 
ing signal is the ellipticity-shear cross orrelation for almost the en- 
tire multipole range, and is typically larger by an order of magni- 
tude compared to the intrinsic ellipticity correlation, too. 

iHevmans et al.l ( l2006t) focus on the quantification of the con- 
tamination of weak lensing measurements due to intrinsic align- 
ments from simulations. They identify haloes in a cosmological n- 
body simulation and model the appearance galaxies in these haloes, 
either by the ellipticity of the halo for an elliptical galaxy or by the 
inclination angle of the angular momentum axis relative to the line 
of sight, allowing for a misalignment between angular momentum 
and symmetry axis of the disk. The derived ellipticity correlation 
function t/^ (r) can be parameterised by a model. 
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Figure 6. Comparison of the weak lensing shear and the intrinsic alignment signal for deep surveys (zmcd = left panel), and for shallow surveys (zmcd =0.1, 
right panel). The intrinsic alignment signal is predicted using a quadratic alignment model with a spin correlation parameter of a = 0.24 (solid line) and 
a = 0.55 (dashed Hne). The long-dashed line indicates the intrins ic lensing signal, with superimposed data point from lensing surveys. The figure was kindly 
provided by R. Crittenden and appears in lCrittenden et a/1 \200l\) . Reproduced with permission from the American Astronomical Society. 
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with two parameters A, B, which agrees well with the alignments 
measurements carried out by Mandelbaum et al. (2006) for the spi- 
rals, but shows deviations in the elliptical galaxy set. Morphologi- 
cally, the intrinsic ellipticity correlation assumes the highest values 
for the elliptical galaxies and much lower values for spirals. 

At the same time, they derive the cosmological weak lens- 
ing signal from the weighted projected surface mass density of the 
simulated density field, and distinguish between intrinsic elliptic- 
ity correlations, which can be dealt with if accurate photometry is 
available by discarding close pairs, and correlations between the 
intrinsic ellipticity and the external gravitational shear. The latter 
correlation can be equally modelled by a simple function, depend- 
ing on the angle of separation 6, 

DjDa, A 



D, e + t 



(100) 



with two parameters A and 6q, and scaling with the lensing effi- 
ciency. As illustrated by Fig. |7] the correlation iff, is smaller com- 
pared to the weak lensing signal by a factor of a few in a deep 
survey with sources at unit redshift, while ellipticals show the 
strongest signal, followed my a realistic mix of morphologies. On 
the contrary, the correlation function is consistent with zero for a 
data set comprising only spiral galaxies. 

The impact of intrinsic alignments on higher- order s t atistic s 
of the lensing signal was investigated by Sembolo ni et al.l 
Using the same alignment model as Hevmans e t al.i ^20061) , they re- 
port a surprisingly strong contamination of the weak shear 3-point 
correlation function. Specifically, they consider the skewness of the 
aperture-mass statistic, where the lensing signal is weighted inside 
a circular aperture, and compare the three possible 3-point func- 
tions between the intrinsic alignment signal and the cosmological 
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Figure 7. Lensing-intrinsic shape cross correlation functions derived from 
;i-body simulations, for spiral galaxies (filled dots), a mixed population 
(triangles) and ellipticals (circles), in comparison to the lensing signal 
(solid line), which is modeled by placing the source galaxies at a redshift 
of unity. The figure w as kindly provided by C. Heymans and appears in 
\Hevmans et al\ \200(jl }. Reproduced with permission from Blackwell Pub- 
lishing. 
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weak lensing signal, based on nummerical data. In shallow surveys 
with median redshifts of Zmd - 0.4, one expects the 3-point corre- 
lation of the intrinsic alignment to be larger than the weak lensing 
by almost an order of magnitude, whereas in medium-deep sur- 
veys wi th Zmci = 0.7 the contam ination amounts to ~ 15%. Conse- 
quently, [Se^^olonrern] ([200^ emphasise the usability of 3-point 
stati stics to test intri nsic alignment models. 

iFu et al. I lioos") estimate the magnitude of the intrinsic align- 
ment contamination on the CFHTLS weak lensing survey, by using 
the m odel for ^^^{0) (see eqn. 1 1001 introduced in iHevmans et al.l 
I2OO6) . Setting the angular scale to do = 1 arcmin, a combined fit 
to the cosmological and intrinsic alignment signal yields the value 
A = 2.2^;|g for the amplitude A, i.e. compatible with zero. In the 
range of angular scales considered, the weak lensing signal domi- 
nates over the intrinsic alignment signal and the contamination is 
weak, limited to be less than 32% and larger than -13%. 

8.6 Suppression of intrinsic ellipticity alignments 

The methods for suppressing intrinsic alignment contamination in 
the cosmological weak shear signal fall roughly into five categories: 
With accurate models for the ellipticity correlations, the intrinsic 
alignment signal can be modelled such that one obtains a prediction 
for C^(^) -I- Citeif) + Ceif) for a given cosmology, and it becomes 
possible to fit the entire observational data set. Sadly, the models 
for the computation of €,:{£) and of the cross term C^^if^) are not 
yet reliable enough. Secondly, one can take advantage of particular 
statistical properties of the intrinsic alignment signal, e.g. the ex- 
citation of vortical patterns in the ellipticity field and statistically 
separate the weak shear signal from the intrinsic alignment signal. 
A third method is the truncation of data, i.e. the rejection of galaxy 
pairs with small line of sight separation, at the expense of reduc- 
ing statistics, especially at small separations. A fourth way of deal- 
ing with the problem is the introduction of new lensing weighting 
functions, which null out the intrinsic signal. At last, it is possible 
to treat the intrinsic alignment signal as a systematic, leave it in the 
data and quantify the bias in the cosmological parameter estimation 
which is due to the presence of intrinsic alignment. Most of the pro- 
posed methods focus on eliminating the contribution Ci({), because 
it was only realised bv Hirat a & Seliak ( 2004), that the C^idC) might 
play an important role. Generally, the methods only assume that the 
intrinsic alignment contamination is local, and can be expected to 
treat (f) better than C^i(£), due to the sh orter correlation length. 

For suppressing intrinsic alignments. ICrittenden et al" ||2003) 
put forward a decomposition of the excitations in the ellipcitity field 
into its gradient and a curl components, which is motivated by the 
fact that gravitational lensing can only excite the gradient compo- 
nent and leaves the curl component untouched, as long as one as- 
sumes the validity of the Bom-approximation and the uncorrelated- 
ness of multiple lensing events along the line of sight. Furthermore, 
the lensing signal is linear in the tidal shear, whe reas the intrinsic 
alignments measure the square of the tidal shear. ICrittenden et al.l 
(2002) show that the two contributions to the ellipticity correlations 
can be statistically disentangled, and emphasise the importance of 
measuring the curl components as a tool for controlling systematics 
in the observation. iMackev etai] l l2002h investigate the ratio of the 
gradient and curl power spectra, and find a typical enhancement 
of the gradient power spectrum over the spectrum of the vortical 
excitations of roughly 3.5 on small angular scales. Hence, the con- 
tamination of the weak lensing spectrum is significantly larger than 
one would naively expect from a determination of the amplitude of 
vortical modes in the ellipticity field due to intrinsic alignment. 



iKing&SchneidejJiooi) propose a weighting scheme for di- 
rectly downweighting close pairs of galaxies, where the shape cor- 
relations are likely due to intrinsic alignments and not of cosmo- 
logical origin. This is justified by the fact that tidal torquing models 
predict only short ranged correlations, due to their dependence on 
the squared tidal field. In the computation of the shear correlation 
function they propose to use a redshift weighting Z(zi,Zj), 

Z(z„z,) = l-exp|-^j, (101) 

Az = Zi - Zj, i.e. a constant weighting with a Gaussian dip of width 
(Tj, which down-weights close pairs in redshift z,. a, corresponds 
to the photometic redshift uncertainty, which is typically larger 
than the intrinsic alignment correlation length. iKing & SchneideJ 
(2002) quantify the loss of cosmological information by discard- 
ing galaxy pairs and the residual contamination using a quadratic 
model motivated by conventional tidal torquing. They conclude that 
the alignment contamination can be effectively suppressed, that the 
observed correlation function looses amplitude by the weighting, 
and that the increase in Poisson-noise due to down-weighting can 
be as large as 25% for large cr-, emphasising the importance of 
good photometry. 

Reduction of intrinsic alignment contamination in the con- 
text of the COMBO- 17 survey, a medium-deep photometric sur- 
vey with median redshift Zm ed - 0.6, is the topic of th e papers by 
iHevmans & Heavens! ('2003) and lHevmans et all ( l2004h : They pur- 
sue the strategy of downweighting galaxy pairs at small separation, 
and improve on lKing & Schneider! ( l2002h by deriving the optimal 
pair weight for a spectroscopic survey, which minimises the errors 
in the weak lensing correlation functino due to intrinsic align ment, 
assuming a specific alignment model ( iHeavens et alj |20^; !jing! 
!2002!) . and due to shot noise. They show that the intrinsic align- 
ment contamination can be reduced to the level of pure shot noise 
for a shallow spectroscopic survey, and that the method applies to 
photometric surveys like COMBO-17 as well. In a further step, they 
are able to constrain the alignment model from data, measure the 
amplitude of the intrinsic ellipticity correlation function and quan- 
tify its impact on the measurement of the cosmological parameters 
n„ and o-g from COMBO-17. 

!King & S chneider (2003) demonstrate that in a tomographic 
measurement with accurate photometry the intrinsic alignment con- 
tamination can be disentangled statistically from the lensing signal, 
and extend their study to show that the degeneracy between Cl„, 
and erg, which is inherent to gravitational lensing and which can be 
partially lifted using redshift estimates for the source galaxies, can 
be dealt with even if intrinsic alignments are present. Technically, 
they prepare artificial data consisting of the weak lensing spectrum, 
generated from the CDM spectrum assuming a ACDM cosmology, 
with an intrinsic ali gnment contamination, using the model taken 
from iHeavens et al.! ( !200Q) . and dilute the data by modeling pho- 
tometry errors in the tomography correlation functions. They con- 
tinue by carrying out a;^--fit to the data with a grid of template 
functions for the lensing signal as well as for the intrinsic align- 
ments and show that the two contributions can be separated, due 
to the fact that the intrinsic alignment causes correlations only on 
short separat ions and sho ws a different redshift dependence. In a 
further studv !King| ( l2005l) extends the analysis to deal with cross 
correlations between intrinsic alignment and the weak lensing sig- 
nal in a linear alignment model, with essentially the same success. 

A very el egant solution to the intrinsic alignment problems 
is proposed by !joachimi & Schneided ( !2008!) : They consider a to- 
mographic setup, where the set of lensed background galaxies are 
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divided into a number of redshift bins and design new line of sight 
weighting functions which cancel out the shear signal originating 
from galaxies in the same tomographic bin, by introduction of a 
line of sight weighting function Biix), which replaces the redshift 
distribution of the galaxies Piix) of conventional lensing tomog- 
raphy. Specifically, the line of sight integral for the weak lensing 
convergence ^, replaces /c, . 



1& Jo a 
with the changed lensing efficiency G{x), 



G,(X) 



dx' B.ix') 



X -X 

X' ' 
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which carries out the nulling of the contribution of the large-scale 
structure inside the tomography bin i, by demanding that the co- 
moving distance Xi the weighting function vanishes, Gixi) = 0, 
assuming locality of the intrinsic alignment effect and should there- 
fore be applicable to the suppression of the C;(E(f)-contribution as 
well, despite its longer correlation length. Quite generally, a loss 
of information is associated with the nulling, and consequently 
one would expect a deterioration in the accuracy of cosmologi- 
cal parameters derived from a measurement of the power spectrum 
Cij(f) of the cross-fluctuat ions in the ^,-field. Noticing the anal- 
ogy to variational problems. ljoachimi & SchneideJ ( l20o3) continue 
the derivation of the weighting functions Bj(x) by constructing the 
Fisher matrix of the spectra Cij{C). The Fisher matrix is defined 
as the expectation value of the second derivative of the logarith- 
mic likelihood with respect pairs of cosmological parameters and 
provides a measure of the sensitivity of the lensing spectra for con- 
straining the cosmological parameters. Specifically, the trace of the 
Fisher matrix, which is equal to of the individual errors cr^ 
of the cosmological parameters due to the Cramer-Rao inequal- 
ity, is maximised, with the nulling condition Gixi) = contained 
in the spectra. They are able to show both analytically and numeri- 
cally that the weighting functions Bi(x) with the desired properties 
can be constructed, that the accuracy loss in parameter estimation 
amounts to typically 30%, and that the parameter estimates remain 
unbiased. 

Considering 3d-lensing measurements, iKitching et al] ( |2008|) 
quantified the impact of intrinsic alignments along with other un- 
certainties such as photometry errors on the estimation of cos- 
mological parameters and in particular the dark energy parame- 
ters, using an extended Fisher-matrix formalism. They describe 
the intrinsic al i gnme nts with the parameterisation proposed by 
iHevmans et al.l 1 2006h and show that uncorrected intrinsic align- 
ment contamination has a strong impact on the estimation of the 
dark energy figure of merit FoM = [AvVoAwCzp)]"', degrading it by 
as much as a factor of 2. Additionally, using biased parameters for 
describing the contamination results in large biases in the estima- 
tion of cosmological parameters, but including prior information on 
the systematics is able to recover the accuracy of the measurement 
to a large extent 

Finally, lAmara & Refregiej koot provide a general frame- 
work for quantifying how systematics in lensing measurements in- 
troduce biased values for the estimated cosmological parameters, 
which can be applied to the intrinsic alignment problem. Their so- 
lution to the intrinsic alignment model would be to measure the 
cosmological parameters from the entire weak lensing data set and 
correct for the estimation bias due to intrinsic alignment model. 



9 SUMMARY AND OUTLOOK 

This review article summarises models for galactic angular mo- 
menta and angular momentum correlations in the cosmological 
large-scale structure, and its implications for gravitational lens- 
ing. Starting from perturbative models and tidal torquing, angular 
momentum distributions and correlation functions are derived and 
simplifications and their applicability discussed. In general, linear 
tidal torquing is a very successful theory for explaining the angu- 
lar momentum magnitude and direction, only the nonlinear stages 
become difficult to describe analytically, as well as the influence of 
baryonic dynamics. The many functional forms used show the am- 
biguity in defining an angular momentum correlation function. In- 
teresting future developments might include higher-order statistics 
of the angular momentum field, which is intrinsically non-Gaussian 
even if it develops from Gaussian initial fluctuations, or for non- 
Gaussianities in the initial conditions, either primordial or due to 
translinear structure formation on small scales. Simulations find 
correlations between angular momenta and the strong tidal fields 
of nonlinear structures, illustrating the richness of the cosmologi- 
cal large-scale structure which can not be quantified using 2-point 
functions alone. 

The primary application of angular momentum models are el- 
lipticity correlations, which have been shown to be a serious con- 
taminant of high-precision weak lensing measurements. A plethora 
of methods has been suggested for removing or suppressing in- 
trinsic ellipticity correlations such that the power of weak lens- 
ing surveys for constraining cosmological parameters can be har- 
nessed. Nevertheless one could expect futher development in pre- 
dicting statistical properties of the intrinsic alignments, especially 
in comparing linear and quadratic alignment models, and numer- 
ical studies of how the ellipticity measured from the luminosity 
distribution is related to the angular momentum. The suppression 
and nulling techniques invented to suppress the intrinsic alignment 
can of course be used to achieve an amplification of the alignment 
signal over the weak lensing, and it might be possible to combine 
higher order image distortions such as weak flexions, which should 
be free of intrinsic alignments, with the weak shear signal, enhanc- 
ing the understanding of intrinsic and weak lensing induced ellip- 
ticity correlations. Insight into the alignment of the stellar light 
distribution with respect to the dark matter host structure could 
po ssibly be g ained using strong lensing (an overview is provided 
bv lMevlan et alJ i2006). by measuring the halo's potential well and 
combining with kinematical data. 
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